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Abstract 



The aim of these lectures is to give an introduction to several topics which lie 
at the intersection of string theory, gravity theory and gravity phenomenology. One 
successively reviews: (i) the "membrane" approach to the dissipative dynamics of 
classical black holes, (ii) the current experimental tests of gravity, and their theoret- 
ical interpretation, (iii) some aspects of the string-inspired phenomenology of the 
gravitational sector, and (iv) some possibilities for observing string-related signals 
in cosmology (including a discussion of gravitational wave signals from cosmic 
superstrings). 



1. Introduction 

The common theme of these lectures is gravity, and their aim is to discuss a few 
cases where string theory might have an interesting interplay either with gravity 
theory, or with gravity phenomenology. We shall discuss the following topics: 

• Classical black holes as dissipative branes. The idea here is to review the 
"classic" work on black holes of the seventies which led to the picture of black 
holes as being analog to dissipative branes endowed with finite electrical resis- 
tivity, and finite surface viscosity. In particular, we shall review the derivation 
of the (classical) surface viscosity of black holes, which has recently acquired a 
new (quantum) interest in view of AdS/CFT duality. 

• Hawking radiation from black holes. To complete our classical account of 
irreversible properties of black holes, we shall also give a direct derivation of the 
phenomenon of Hawking radiation, because of its crucial importance in fixing 
the coefficient between the area of the horizon and black hole "entropy". 

• Experimental tests of gravity. Before discussing possible phenomenological 
consequences of string theory in the gravitational sector, we find useful to sum- 
marize the present status of experimental tests of gravity, as well as the theoretical 
frameworks used to interpret them. In particular, we emphasize that binary pulsar 
experiments have already given us accurate tests of some aspects of strong-field 
(and radiative) relativistic gravity. 

• String-inspired phenomenology of the gravitational sector In this section we 
shall discuss (without any attempt at completeness) some of the ideas that have 
been suggested about observable signals possibly cormected to string theory. In 
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particular, we shall discuss the cosmological attractor mechanism which leads 
to a rather rich gravitational phenomenology that will be probed soon by various 
gravitational experiments. 

• String-related signals in cosmology. After discussing a few alternatives to 

slow-roll inflation (and the possible relaxation of the Lyth bound when using 
non-linear kinetic terms for the inflaton), one discusses in some detail cosmic 
superstrings. We explain, in particular, how one computes the gravitational wave 
burst signal emitted by the cusps that periodically form during the dynamical 
evolution of generic string loops. 

A final warning: by lack of time (and energy), no attempt has been made to 
give exhaustive and fair references to original and/or relevant work. The given 
references are indicative, and should be viewed as entry points into the relevant 
literature. With the modern, web-based, easy access to the scientific hterature it 
is hoped that the reader will have no difficulty in using the few given references 
as starting points for an instructive navigation on the vast sea of the physics liter- 
ature. 



2. Classical black holes as dissipative branes 



Early work on (Schwarzschild, Reissner-Nordstrom, or Kerr-Newman) black 
holes (BHs) in the 1950's and 1960's treated them as passive objects, i.e., as 
given geometrical backgrounds (and potential wells). This viewpoint changed in 
the early 1970's when the study of the dynamics of BHs was initiated by Pen- 
rose [1], Christodoulou and Ruffini [2, 3], Hawking [4], and Bardeen, Carter and 
Hawking [5]. In the works [1-5], only the global dynamics of BHs was consid- 
ered, i.e.. their total mass, their total angular momentum, their total irreducible 
mass, and the variation of these quantities. This viewpoint further evolved in the 
works of Hartle and Hawking [6], Hanni and Ruffini [7], Damour [8-10], and 
Znajek [1 1 1, which studied the local dynamics ofBH horizons. In this new ap- 
proach (which was later called the "membrane paradigm" [12]) a BH horizon is 
interpreted as a brane with dissipative properties, such as, for instance, an elec- 
trical resistivity p, equal to 377 Ohms [8, 11] independently of the type of BH, 
and a surface (shear) viscosity, equal to r; = [9, 10]. When divided by the 
entropy density found by Hawking (S/A = j), the latter shear viscosity yields 
the ratio j^, a result which has recently raised a renewed interest in connection 
wifli AdS/CFT, flirough the work of Kovtun, Son, and Starinets [13, 14]. 
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2.1. Global properties of black holes 

Let us start by reviewing the study of the global dynamics of BHs. Initially, BHs 
were thought of as given geometrical backgrounds. In the case of a spherically 
symmetric object of mass M without any additional attribute, Schwarzschild de- 
rived the first exact solution of Einstein's equations only a few weeks after Ein- 
stein had obtained the final form of the field equations. Schwarzschild's solution 
is as follows. In 3 + 1 dimensions, setting G — c — \, the metric for a spherically 
symmetric background can be written in the form 

ds^ = -A{r)dT'^ + B{r)dr'^ + {dO"^ + sm^OdLp^) (2.1) 

where T denotes the usual Schwarzschild-type time coordinate, and where the 
coefficients A (r) and B (r) read 

Mr) ^ ^^^^^ 



B{r) 



Mr) 



This result was generalized in independent works by Reissner, and by Nordstrom 
(1918) for electrically charged spherically symmetric objects, in which case A (r) 
and B (r) are given by 



A{r) = 1-^ + ^ 
B{r) - 



(2.3) 



Mr) 



We shall not review here the long historical path which led to interpreting the 
above solutions, as well as their later generalizations due to Kerr (who added to 
the mass M the spin ^Ql, and Newman et al. (mass, spin and charge), as BHs. Up 
to the 1960's BHs were viewed only as passive gravitational wells. For instance, 
one could think of adiabatically lowering a small mass m at the end of a string 
until it disappears within the BH, thereby converting its mass-energy mc^ into 
work. More realistically, one was thinking of matter orbiting a BH and radiating 
away its potential energy (up to a maximum, given by the binding energy of the 
last stable circular orbit around a BH). This viewpoint changed in the 1970's, 
when BHs started being considered as dynamical objects, able to exchange mass, 
angular momentum and charge with the external world. Whereas in the simplest 
case above, one uses the attractive potential well created by the mass M without 
extracting energy from the BH, Penrose [1], showed that energy could in princi- 
ple be extracted from a BH itself by means of what is now called a (gedanken) 
"Penrose process" (see FIG. [T]i. Namely, if one considers a time-independent 

'Note that in the case of a spinning BH, one often introduces the useful quantity a = J/M, i.e., 
the ratio of the total angular momentum to the mass of the BH, which has the dimension of length. 
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Figure 1. In this figure, we schematically illustrate the "Penrose process", i.e., the splitting of an 
ingoing particle into one that falls into the BH and another that exits at infinity. 

background and a BH that is more complicated than Schwarzschild's, say a Kerr 
BH, one may extract energy using a test particle 1 coming in from infinity with 
energy Ex, angular momentum p^^, and electric charge ei. By Noether's theo- 
rem, the time-translation, axial and U{\) gauge symmetries of the background 
guarantee the conservation of E, and e during the "fall" of the test particle. 
Moreover, if, in a quantum process, the test particle 1 splits, near the BH, into 
two particles 2 and 3, with E^, Pip^, 62, and £3, p^p^, 63 respectively, then, under 
certain conditions, one finds that particle 3 can be absorbed by the BH, and that 
particle 2 may come out at infinity with more energy than the incoming particle 
1. A detailed analysis of the efficiency of such gedanken Penrose processes by 
Christodoulou and Ruffini [2, 3] then led to the understanding of the existence of 
a fundamental irreversibility in BH dynamics, and to the discovery of the BH 
mass formula. Let us explain these results. 

The basic idea is to explore the physics of BHs through a sequence of infinites- 
imal changes of their state obtained by injecting in them some test particles. One 
starts by writing that the total mass-energy, spin and charge of the BH change, by 
absorption of particle 3, as 



6M 
6J 
6Q 



Es = El — E2, 
J3 = "^1 — J2, 
63 = ei - 62. 



(2.4) 
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This preliminary result can be further exploited by making use of the Hamilton- 
Jacobi equation. Considering an on-shell particle of mass /i, and adopting the 
( — h ++) signature, the Hamilton-Jacobi equation reads 

g^''(p^-eA^)(p,-eA,) - V, (2.5) 

in which = dS/dx^, S is the action and the partial derivatives are taken 
w.r.t. the coordinates x^. The details of the splitting process will be irrelevant, 
as only particle 3 matters in the calculation. In an axially symmetric and time- 
independent background, S can be taken as a linear function of T and (p, 

S^-ET + p^if + S{r,d). (2.6) 

where E — —px — —po is the conserved energy, p^p is the conserved (p- 
component of angular momentum and the last term is the contribution from terms 
that depend on the angle 9 and on the radial distance r. Let us consider the case 
of a Reissner-Nordstrom BH, where calculations are easier: the inverse metric is 
easily computed and (|2.5[) can then be written explicitely as 



^ (po - eAof + A{r)pl + \ ( p^ + -^pl ] = -fi^ (2.7) 



A{r) \ sm^ 

which we re-write as 

/ r2 

\2 _ /l/'„^2„2 I I ,,2 



{po - eAor = A{ryp'^ + A{r) ^fi' + — j (2.8) 

The electric potential is —Aq = +V = +Q/r. The above expression is quadratic 
in E (it is the generalization of the famous flat-spacetime E^ = /^^ + p^) and one 
finds two possible solutions for the energy as a function of momenta and charge 
(seeFIG.|51i: 




E=^±\lA{r)^p^ + A{r)(^fi^ + ^y (2.9) 

In flat space, A (r) — 1, so that, if we ignore charge, we recover the usual Dirac 
dichotomy on the choice of the + or — sign between particle and antiparticle: 
E = -^\/ IJ? + p^. This shows that one should take the plus sign in the equation 
above. We remind the reader that for a charged BH, there exists a regular horizon 
only if Q < M (which can be interpreted as a BPS bound). [We have set G = 1]. 
Remembering that A{r) = 1 — 2M/r + Q"^ jr^, there exists both an outer and 
an inner horizon defined by r± = M ± \J A'P — (which are the two roots 
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Figure 2. This figure depicts the classically allowed energy levels (shaded region) as a function 
of radius, for test particles in the neighborhood of a BH. There exist positive- and negative-energy 
solutions, corresponding (after second quantization) to particles and anti-particles. Classically (as in 
the Penrose process) one should consider only the "positive-square-root" energy levels, located in the 
upper shaded region. The white region is classically forbidden. Note the possibility of tunneling (this 
corresponds to particle creation via the "super-radiant", non-thermal mechanism briefly mentioned 
below). 



of A{r) = 0). The horizon of relevance for BH physics is the outer one r+ = 
M + ^/W^^^^Cp (it gives the usual result 2M when Q = 0). As particle 3 
is absorbed by the BH, we can compute its (conserved) energy when it crosses 
the horizon, i.e., in the Umit where the radial coordinate r is equal to r+. This 
simplifies the expression of to 

= — + b1, (2.10) 
r+ 

where we have introduced the contravariant component = cf^Pr = A{r)pr, 
which has a finite hmit on the horizon. Note the presence of the absolute value 
of p'^ (coming from the limit of a positive square-root). The change in the mass of 
the BH is equal to the energy Es of the particle absorbed, i.e., particle 3. Using 
^3 = SQ, this yields 



From the positivity of \p'^\ we deduce that 



(2.12) 
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We have derived an inequality and have thereby demonstrated (by following 
Christodoulou and Ruffini) the UTeversibility property of BH energetics. There 
exist two types of processes, the reversible ones with an '=' sign in ( 2.12| i, 



and the irreversible ones with an '>' sign. The former ones are reversible be- 
cause if a BH first absorbs a particle of charge +e with vanishing \p^\ (thereby 
changing its mass by 6'M = eQ/r^{M,Q) and its charge by 5'Q = e), and 
then a particle of charge — e with vanishing \p^\ (thereby changing its mass by 
S"M = -eQ/r+{M,Q) and its charge by S"Q = -e), it will be left, at the 
end, in the same state as the original one (with mass M + S'M + 5"M — M 
and charge Q + 5'Q + 5"Q = Q). Evidently, such reversible transformations 
are delicate to perform, and one expects that irreversibility will occur in most BH 
processes. The situation here is clearly similar to the relation between reversible 
and irreversible processes in thermodynamics. 

The same computation as for the Reissner-Nordstrom BH can be performed 
for the Kerr-Newman BH. One obtains in that case, by a slightly more compli- 
cated calculation, 

a5J + r+Q5Q + aWg 

5M T^--2 = 2,2 — \P ■ (2-13^ 

in which r+(M, J, Q) = M + - - a"^. We recall that a = J/M, and 

that one has the bound + (J/M)^ < M^. 

The idea now is to consider an infinite sequence of infinitesimal reversible 
changes (i.e., 0), and to study the BH states which are reversibly connected 
to some initial BH state with given mass M, angular momentum J and charge 
Q. This leads to a partial differential equation for 6M, 

^^j^a_5J+r^^ (2.14) 

which is found to be integrable. Integrating it, one finds the Chris todoulou- 
Ruffini mass formula [3] 



Here the irreducible mass Mjrr = 5 y ' + 
stant. The mass squared thus appears as a function of three contributions, with 
one term containing the square of the sum of the irreducible mass and of the 
Coulomb energy, and the other one containing the rotational energy. Inserting 



this expression into Eq. (2.13 1, one finds 



6Mi„ > 



(2.16) 
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with JA/ii-i = under reversible transformations and JMin > under irre- 
versible transformations. The irreducible mass A/jn can only increase or stay 
constant. This behaviour is certainly reminiscent of the second law of thermody- 
namics. The free energy of a BH is therefore M — Mi„, i.e., this is the maximum 
extractable energy. In this view, BHs are no longer passive geometrical back- 
grounds but contain stored energy that can be extracted. Actually, the stored en- 
ergy can be enormous because a BH can store up to 29 % of its mass as rotational 
energy, and up to 50 % as Coulomb energy! 

The UTeducible mass is related to the area of the horizon of the BH, by ^ = 
IBttM^j. so that in a reversible process 5 A — 0, while in an irreversible one 
5 A > 0. Hawking showed [4] that this irreversible evolution of the area of the 
horizon was a general consequence of Einstein's equations, when assuming the 
weak energy condition. He also showed that in the merging of two BHs of area 
Ai and A2, the total final area satisfied Atot > ^1 + ^2- 

Such results evidently evoque the second law of thermodynamics. The analog 
of the first law [di?(S', extensive parameters) = dW + dQ, where the work 
dW is linked to the variation of extensive parameters (volume, etc.) and where 
dQ = TdS is the heat exchange] reads, for BH processes, 

dM{Q,J,A)^VdQ + ndJ+—dA. (2.17) 

Comparing this result with expression ( |2.13| l, one has 



V 

n 



a 



and 

which, in the Kerr-Newman case, is given by 



(2.18) 



(2.19) 



(2.20) 



V is interpreted as the electric potential of the BH, and Q, as its angular velocity. 



Expression (2.17 1 resembles the usual form of the first law of thermodynamics in 
which the area term has to be interpreted as some kind of entropy. The parameter 
g is called the "surface gravity". [In the Schwarschild case, it reduces to AI /r\ 
(in G = 1 units), i.e., the usual formula for the surface gravitational acceleration 
g — GM/R^.]. In the Les Houches Summer School of 1972, a more general 



String theory, gravity and experiment 



11 



version of the first law was derived, that included the presence of matter around 
the BH, and energy exchange [5]. An analog of the zeroth law was also derived 
[15], in the sense that the surface gravity g (which is analog to the temperature) 
was found to be uniform on the surface of a BH in equilibrium. 

In 1974, Bekenstein went further in taking seriously (and no longer as a simple 
analogy) the thermodynamics of BHs. First, note that one can write the formal 
BH "heat exchange" term in various ways 

dQ = TdS = ^ = 4gMi„dMi„. (2.21) 
Stt 

In light of this, is the appropriate physical analog of the entropy the irreducible 
mass or the area of a BH? Is the analog of temperature proportional to the sur- 
face gravity g or to the product Mjrr gl Can one give a physical meaning to the 
temperature and entropy of a BH ? To address such questions, Bekenstein used 
several different approaches. 

In particular, he used Carnot-cycle-type arguments. For instance, one may ex- 
tract work from a BH by slowly lowering into it a box of radiation of infinitesimal 
size. In fact, in this ideal case, one can theoretically convert all the energy of the 
box of radiation, mc^, into work. The efficiency of Camot cycles is defined in 
terms of both a hot and a cold source as 

^7=1-^. (2.22) 

-ihot 

From what we just said, it would seem that the efficiency of classical BHs as 
thermodynamic engines is 100%, rj = \. This would then correspond to a BH 
temperature (= the cold source) Tbh = Tcoid = 0- The point made by Beken- 
stein was that this classical result will be modified by quantum effects. Indeed, 
one expects (because of the uncertainty principle) that a box of thermal radiation 
at temperature T (made of typical wavelengths A ~ 1/T) cannot be made in- 
finitesimally small, but will have a minimum finite size ~ A. From this Umit on 
the size of the box, Bekenstein then deduced an upper bound on the efficiency ry, 
and therefore a lower bound on the BH temperature Tbh 0. 

Let us indicate another reasoning (of Bekenstein) which suggests that the ab- 
sorption of a single particle by a BH augments its surface by a finite amount 
proportional to h. As we said above the change of BH energy as it absorbs a 
particle is (when a = 0, for simpUcity) 

£;3 = ^+lim^^^^|pn (2.23) 
r+ 

We also showed that the transformation will be reversible (i.e., will not increase 
the surface area of the BH) only if limr^r+ \ = 0. However, for this to be 
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trae both the (radial) position and the (radial) momentum of the particle must be 

exactly fixed: namely, r = r+ and = 0. This would clearly be in contradiction 
with the Heisenberg uncertainty principle. Technically, we must consider the 
conjugate momentum to the position r which is the covariant component Pr of 
the radial momentum (instead of the contravariant component p^ used in the 
equation above). The uncertainty relation therefore reads 

5r5pr >]^h. (2.24) 

Near the horizon (i.e. when 5r = r — r+ is small), the contravariant radial 
momentum reads (using g^'^ = l/g^^ = A{r)) 

p'' = A{r)pr 

(rj-r+HTj-r_) 

Srpr, 



1"+ 



so that Heisenberg's uncertainty relation yields a lower bound for p^. We can 

reexpress this lower bound in terms of the BH surface gravity g introduced above 
by noting that the partial derivative of A w.r.t. r, (f^)^^, entering the last equa- 
tion, is proportional to g: 

'dA\ 



2g. (2.26) 
This then gives 

p'' ~ 2g6r6pr > gh (2.27) 
>From the relation 5M = ^^'^m^q) + \p''\r+, we finally obtain 

SM-^ = \p^\>gh, (2.28) 
r+ 

which can be rewritten as 

6A > Snh. (2.29) 

In other words, quantum mechanics tells us that when one lets a particle fall into 
a BH, one cannot do so in a perfectly reversible way. The area must increase 
by a quantity of order h. If (still following Bekenstein) one considers that the 
irreversible absorption of a particle by a BH corresponds to the loss of one bit 
of information (for the outside world), we are led to the idea of attributing to 
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a BH an entropy (in the sense of "negentropy") equal (after re-introducing the 
constants c and G) to [16] 

Sbh = a— A, (2.30) 

hG 

with a dimensionless numerical coefficient equal to a = ln2/87r according to 
the reasoning just made. More generally, Bekenstein suggested that the above 
formula should hold with a dimensionless coefficient a « O (1), without being 
able to fix in a unique, and convincing, manner the value of a. This result in turn 
imphes (by applying the law of thermodynamics) that one should attribute to a 
BH a temperature equal to 

TBH = ^-g. (2.31) 
8na c 

This attribution of a finite temperature to a BH looked rather strange in view of 
the definition of a BH has being "black", i.e., as allowing no radiation to come out 
of it. In particular, Stephen Hawking resisted this idea, and tried to prove it wrong 
by studying quantum field theory in a BH background. However, much to his own 
surprise, he so discovered (in 1974) the phenomenon of quantum radiation from 
BH horizons (see below) which remarkably vindicated the physical correctness 
of Bekenstein's suggestion. Hawking's calculation also unambiguously fixed the 
numerical value of d to be d = j [17]. [We shall give below a simple derivation 
(fromRef. [18]) of Hawking's radiation.] 

Summarizing so far: The results on BH dynamics and thermodynamics of the 
early 1970's modified the early view of BHs as passive potential wells by endow- 
ing them with global dynamical and thermodynamical quantities, such as mass, 
charge, irreducible mass, entropy, and temperature. In the following section, we 
shall review the further changes in viewpoint brought by work in the mid and 
late 1970's ( [6, 8-11]) which attributed local dynamical and thermodynamical 
quantities to BHs, and led to considering BH horizons as some kind of dissi- 
pative branes. Note that, in the following section, we shall no longer consider 
only Kerr-Newman BHs (i.e., stationary BHs in equilibrium, which are not dis- 
torted by sources at infinity). We shall consider more general non-stationary BHs 
distorted by outside forces. 

2.2. Black hole electrodynamics 

The description of BHs we give from here on is essentially "holographic" in 
nature since it will consist of excising the interior of a BH, and replacing the 
description of the interior BH physics by quantities and phenomena taking place 
entirely on the "surface of the BH" (i.e., the horizon). The surface of the BH 
is defined as being a null hypersurface, i.e., a surface everywhere tangent to the 
lightcone, separating the region inside the BH from the region outside. As just 
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said, we ignore the region inside, including the spacetime singularity, and con- 
sider the physics in the outside region, completing it with suitable "boundary 
effects" on the horizon. These boundary effects are fictitious, and do not really 
exist on the BH surface but play the role of representing, in a holographic sense, 
the physics that goes on inside. In the end, we shall have a horizon, a set of 
surface quantities on the horizon and a set of bulk properties outside the horizon. 
We first consider Maxwell's equations, namely F^^ = d^Ai, — d^A^, and 

W.F^^" = 47^J^ 

v^j^ = 0. ^^-^^^ 

A priori, the electromagnetic field Fj^i, permeates the full space time, existing 
both inside and outside the horizon, and the current, i.e., the source term of 
that carries charge, is also distributed both outside and inside the BH. In order to 
replace the internal electrodynamics of the BH by surface effects, we replace the 
real Ff^i^{x) by _F^„^(.t)Gh, where Qh is a Heaviside-like step function, equal to 
1 outside the BH and inside. Then we consider what equations are satisfied by 
this Gjj-modified electromagnetic field. The corresponding modified Maxwell 
equations contain two types of source terms, 

47r(J^e + j^^), ^^-'^^ 
where we have introduced a BH surface current j'^ as 

= ^^""V.e. (2.34) 

This surface current contains a Dirac (5-function which restricts it to the horizon. 
Indeed, let us consider a scalar function (x) such that ^p{x) =Q on the horizon, 
with (f {x) < inside the BH, and (f {x) > outside it. The BH 0-function in- 
troduced above is simply equal to Oh = S{f (a;)), where denotes the standard 
step function of one real variable. Therefore, the gradient of 6h reads 

d^Qn = d,,e {x)) = 6 (ip (x)) d^,p, (2.35) 

where 6 is the (one dimensional) usual Dirac delta, so that 6 {(f (x)) is a delta 
function with support on the horizon. Morally, the gradient yields a vector 
"normal to the horizon". In the case of a BH (by contrast to the usual case of 
a hypersurface in Euclidean space), there exists an extra subtlety in the exact 
definition of the normal to the horizon. The horizon is a null hypersurface which 
by definition is normal to a null covariant vector satisfying both ifj.'^ = and 
l^dx^ for any infinitesimal displacement dx'' within the hypersurface. Since 
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is null, it cannot be normalized in the same way as in Euclidean space. This 
leads to an ambiguity in the physical observables related to i^^. In stationary- 
axisymmetric spacetimes, one uniquely normalizes £^ by demanding that the 
corresponding directional gradient I'^df^ be of the form d/dt + i^d/dcj) (with 
a coefficient one in front of the time-derivative term). We shall assume (in the 
general non-stationary case) that is normalized so that its normalization is 
compatible with the usual normalization when considering the limiting case of 
stationary-axisymmetric spacetimes. Anyway, given any normalization, there 
exists a scalar uj such that 

e^j, = ujd^ip, (2.36) 
and we can then define an "horizon (5-function" 

5h=-5{^), (2.37) 

UJ 

such that 

d^Qn = it^Sn- (2.38) 
This leads to defining a "BH surface current density" 

K^' = —F'^'e^. (2.39) 
An 

With this definition, the BH current reads 

= X^fe, (2.40) 

and satisfies 

W^{eHJ^ + Kf'SH)^0, (2.41) 

which is a conservation law for the sum of the outside bulk current QrJ^ and 
of the boundary current K^Sh- In picturesque terms, the surface current K'^Sh 
effectively "closes" the external current lines penetrating the BH (analogously to 
the case of external currents being injected in a perfect conductor and leading to 



currents flowing on its surface). In addition, Eq. (2.39 1 shows that this surface 
current is linked to the electromagnetic fields which are on the horizon. We have 
thus endowed the horizon with surface quantities, defined uniquely and locally 
on the horizon. 

Before we proceed, we introduce a convenient coordinate system to describe 
the physics on the horizon of a general BH. We assume some regular "slicing" of 
the horizon and its neighbourhood by some (advanced) Eddington-Finkelstein- 
like time coordinate t ~ x'^. Then we assume that the first coordinate is 
such that it is equal to zero on the horizon (like r — r+ in the Kerr-Newman 
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case). Finally for A= 2,3 denote some angular-like coordinates on the two- 
dimensional spatial slice St = t) of the horizon. In this coordinate system, 
we normalize such that 

^^^. = 1+^^^- (2.42) 

Here, we have used the fact that the "normal" vector i^, being null, is also tan- 
gent to the horizon, so that ^^9^ is a general combination of d/dt and d/dx^ 
but has no component along the "radial" (or "transverse") coordinate x^. Be- 
cause i^^ is a vector tangent to the hypersurface, we can consider its integral lines 

= dx^ /dt, which lie within the horizon. These integral curves are called the 
generators of the horizon. They are null geodesies curves, lying entirely within 
the horizon. 

Expression i2A2i for the directional gradient along £^ suggests that be 



interpreted as the velocity of some "fluid particles" on the horizon, which are 
the "constituents" of a null membrane. Similarly to the usual description of the 
motion of a fluid, one has to keep track of the changes in the distance between two 
fluid particles as the fluid expands and shears. For a usual fluid, one considers the 
gradient of the velocity field, splitting it into its symmetric and anti-symmetric 
parts. The antisymmetric part is simply a local rotation which has no incidence 
on the physics and can be ignored. The symmetric part is further split into its 
trace and tracefree parts, namely 

^ {d,Vj + djVi) = (Jij + ■ v6ij (2.43) 

where d is the spatial dimension of the considered fluid (which will be d = 2 in 
our case). Here the first term describes the shear, and the second describes the 
rate of expansion. We will see later how the BH analogs of these quantities are 
defined. For the moment let us consider the distances on the horizon. They are 
measured by considering the restriction to the horizon of the spacetime metric 
(which is assumed to satisfy Einstein's equations). As we are considering a null 
hypersurface, we have 

ds^Ui^o = lAB {t,x'^) {dx^ - v'^dt) {dx^ - v^dt) (2.44) 

where v"^ = ^g^. Note that ds^ is a degenerate metric: indeed, on a (three- 
dimensional) null hypersurface, there is no real time direction (ds^ vanishes 
along the generators). One has only two positive-definite space dimensions along, 
e.g., the spatial slices St- This metric describes the geometry on the horizon from 
which one can compute the area element of the spatial sections St 

dA = v/dctTAsdx^ A da;^. (2.45) 
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One can decompose the current density K'^ into a time component an = K^, 
and two spatial components tangent to the spatial shces St (t = const.) of 
the horizon, 

K^'d^ = aHdt + K^dA (2.46) 
in which dt = t^d^ — v^Oa so that 

K^d^ = ani" + {K^ - (7HV^)dA (2.47) 

The total electric charge of the spacetime is defined by a surface integral at oo, 
say 

Qtot = ^^^ Ipf^^dS^,. (2.48) 

This result can be re-written as the sum of a surface integral on the horizon and 
a volume integral in between the horizon and oo. The volume integral is simply 
the usual charge contained in space, so that we can define the BH charge Qh as 

QH = ^j J^F^^dS^,, (2.49) 

where the tensorial horizon surface element reads dS'^,y = ^e^,ypcrdx'' A Ax" = 
{n-ix^v — n^^n) dA. Here, is a second null vector, which is transverse to the 
horizon, and which is orthogonal to the spatial sections St. It is normahzed such 
that n^^^ = +1. Using the definitions above of the BH surface current, one 
easily finds that the total BH charge can be rewritten as 

Qh = f^'^HdA, (2.50) 

where aH is the time component of the BH surface current introduced above. 

Though it is a priori only the integrated BH charge which has a clear physical 
meaning, it is natural to consider the density an appearing in the above surface 
integral as defining a charge distribution on the horizon. Then the link 

1 

can be thought of as being analog to the result a = -^E^rii giving the electric 
charge distribution on a metallic object. This can again be viewed as part of a 
holographic approach in which the interior of the BH is replaced by boundary 
effects. This analogy extends to the (spatial) currents flowing along the surface 
of the BH. Indeed, using the conservation law (Ojj J** + K^^Sh) = 0, which 
is just a Bianchi identity, one has 



an = K^'n^ = —F^'^n^l, (2.51) 
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This shows, in a mathematically precise way, how an external current injected 
"normally" to the horizon "closes" onto a combination of currents flowing along 
the horizon, and/or of an increase in the local horizon charge density. One can 
also introduce the electromagnetic 2-form and restrict it on the horizon. It then 
defines the electric and magnetic fields on the horizon according to 

^Ff.^dx^' A dx^H = EAdx"^ A d< + Bj_dA. (2.53) 

Taking the exterior derivative of the left-hand-side then gives 

V xE=-^dti^jB^). (2.54) 

which relates the electric and magnetic fields on the horizon. 

>From the various formal definitions above, one also gets the following rela- 
tion 

Ea + eAsB^v^ - 4^7As {K^ - <Jhv^) , (2.55) 

or 

E + V X Bj_ = in (^K - auv^ ■ (2.56) 

We recognize here a BH analog of the usual Ohm's law relating the electric field 
to the current (especially in the case where v 0, i.e., in the absence of the 
various "convection effects" linked to the horizon "velocity" v). From this form 
of Ohm's law, we can read off that BHs have a surface electric resistivity equal 
tO/9 = 47r = 377 Ohm [8, 11]. 

Let us give an example in which this BH Ohm's law can be "applied" to a 
specific system. We consider for simplicity the case of a Schwarzschild BH and 
set up an electric circuit "on the surface of the BH" by injecting on the North pole 
(through an electrode penetrating the horizon under a polar angle 6i, with, say, 
6*1 <C 1) an electric current /, and letting it escap^from the South pole (via an 
electrode penetrating the horizon under a polar angle 62, with, say, tt — 62 ^ 1). 
When viewing the BH as a membrane with surface resistivity p, this set up will 
give rise to a fictitious current flow on the horizon, closing the circuit between 
the North and the South poles. Associated to the current flow on the horizon, 
there will be a potential drop V between the poles. This potential drop is simply 
given by the usual Ohm's law, V = RI, i.e., the product of the current / by a 
"resistance" R: 

V = -Aoiei) + Ao{92) = RI. (2.57) 



-Actually, as (classical) charges cannot escape from a BH, we need to inject in the South electrode 
a flow of negative charges (while injecting a flow of positive charges down the North pole). 



String theory, gravity and experiment 



19 



The BH resistance R can be computed in two different ways, either by solving 
Maxwell's equations in a Schwarzschild background, or by computing, in usual 
Euclidean space, the total resistance of a spherical metallic shell with a uniform 
surface resistivity p = An (by decomposing the problem in many elementary 
resistances, some being in parallel, and others in series). Both methods give the 
same answer, namely 

tan^ 

i? = 21n 1-, (2.58) 

tan^ 

expressed in units of 30 fi. [^This result is saying that the typical total resistivity 
of a BH is of the order of 30 il. In addition, if one considers a rotating BH placed 
in a magnetic field out of alignment with its axis of rotation (a field uniform at oo, 
but distorted on the horizon), one expects to find eddy currents on the horizon, 
currents which dissipate the energy. These currents exist, can be computed and 
do indeed brake the rotation of the BH. In such a situation, one also finds a torque 
which acts to restore the alignment of the BH with the field [8]. 

2.3. Black hole viscosity 

In the previous section, we introduced the electromagnetic dissipative properties 
of a BH, using a holographic approach which kept the physics outside up to in- 
finity, and replaced the physics inside the BH by defining suitable quantities on 
the horizon, and then showed that they satisfied equations similar to well-known 
ones (such as Ohm's law). We now turn to the viscous properties of BHs and 
show how suitably defined "surface hydrodynamical" quantities satisfy a sort of 
Navier-Stokes equation. Technically, we would like to do, for the gravitational 
surface properties, something similar to what we did for the electrodynamic prop- 
erties. Namely, we would like to replace the spacetime connection, say uj, by 
some sort of "screened connection" Qh^^, and see what kind of quantities and 
physics will be so induced on the surface of the BH. However, Einstein's equa- 
tions being nonlinear, one cannot simply use a BH step function Qh as was done 
for BH electrodynamics. We shall therefore motivate the definition of suitable 
"surface quantities" related to in a slightly different way and then study the 
evolution of these surface quantities and their connection to the physics outside 
the horizon, up to cx). Our presentation will be sketchy; for technical details, 
see [9,10,19,20]. 

Let us start by considering an axisymmetric spacetime. Then there exists a 
Killing vector m = m^^d/dx^ = d/d(p, to which, by Noether's theorem, one 

^Indeed, in CGS-Gaussian units (as used, say, in the treatise of Landau and Lifshitz) 30 ohms is 
equal to the velocity of light (or its inverse, depending on whether one uses esu or emu). Then, when 
using (as we do here) units where c = 1, 30Q = 1. 
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can associate a conserved total angular momentum, which can be written as a 
surface integral at oo. The total angular momentum w.rt. y reads 



where dS^v = ^s^vpadxP A da;'^, V denotes a covariant derivative, and the 
surface integral is performed over the 2-sphere, 5oo. This starting point is the 
analog of the surface-integral expression for the total electric charge used above 
to motivate the definition of a BH surface charge distribution. 

In a way similar to what was done in the electromagnetic case, we can use 
Gauss'theorem to rewrite this integral as the sum of two contributions: (i) a vol- 
ume integral (over the 3-volume contained between the horizon and infinity) mea- 
suring the angular momentum of the matter present outside the horizon, and (ii) 
a surface integral over a (topological) 2-sphere Sh at the horizon, representing 
what we can call the BH angular momentum Jh, i.e., 



where Jh is given by the same surface-integral formula as Jqo, except for the 
replacement of S'oo by Sh as integration domain. 

The horizon being tangent to the lightcone, one defines on the horizon, as 
above, a null vector both normal and tangent to it. can in turn be comple- 
mented by another null vector such that ^'^n^ = 1 and such that the surface 
element dS^i, can then be re-expressed as {n^£^ — n,y£f^) dA. Remembering that 
the Killing symmetry preserves the generators of the horizon i.e., the commuta- 
tor [£, m] = 0, one has P'Vvm'^ = myS/^l^, so that we can re-express the BH 
angular momentum Jh as the following surface integral 



This result involves the directional (covariant) derivative of the horizon normal 
vector I along a vector m which is tangent to the horizon. The crucial point now 
is to realize that, very generally, given any hypersurface, the parallel transport 
along some tangent direction, say t, of the (normalized) vector i normal to the 
hypersurface yields another tangent vector. The technical proof of this fact con- 
sists of starting from the fact that I ■ I = e, where e is a constant which is equal 
to ±1 in the case of a time-like or spacelike hypersurface, and to in the case (of 
interest here) of a null hypersurface. Then, taking the directional gradient of this 
starting equality along an arbitrary tangent vector t yields (V(€).^ = 0. From 
this result, one deduces that the vector (Vj-£) must be tangent to the hypersur- 
face. Therefore, there exists a certain linear map K, acting in the tangent plane 




(2.59) 




(2.60) 




(2.61) 
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to the hypersurface, such that Vp^ = K{T). For a usual (time-Uke or space-like) 
hypersurface, the linear map K is called the "Weingarten map" and is simply the 
mixed-component -fCj version of the extrinsic curvature of the hypersurface (usu- 
ally thought of a being a symmetric covariant tensor Kij). On the other hand, in 
the case of a null hypersurface, there is no unique way to define the analog of 
the covariant tensor K^j (where the indices i,j are "tangent" to the hypersur- 
face), but it is natural, and useful, to consider the mixed-component tensorifj, 
intrinsically defined as the Weingarten map K in V^^ = K{^). 

To explicitly write out the various components of the hnear map K (acting 
on the hypersurface tangent plane), we need to define a basis of vectors tangent 
to the horizon. This basis contains the null vector I (which is both normal and 
tangent to the horizon), and two spacelike vectors. Using a coordinate system 
{A = 2,3) of the type already introduced (with the horizon being 
located at = 0), we can choose, as two spacelike horizon tangent vectors, the 
vectors ba = Oa- Then one finds that the Weingarten map K is fully described 
by the set of equations 



(2.62) 



The first equation follows from the fact that I is tangent to a null geodesic ly- 
ing within the null hypersurface. [In turn, this follows from the fact that ^ is 
proportional to the gradient of some scalar, say (satisfying the eikonal equa- 
tion (Viys)^ = 0).] The coefficient g entering the first equation defines (in the 
most general manner) the surface gravity of the BH. We see that it represents 
one component of the Weingarten map K. The other components are the two- 
vector Qai and the mixed two-tensor D^. One can show that the components 
D'^ are the mixed components of a symmetric two-tensor Dab, which mea- 
sures the "deformation", in time, of the geometry of the horizon. We remind 
the reader of the expression of the horizon metric, introduced above, d.s^jif = 
7ab {t, x) {dx^^ — v^dt) {dx^ — v^dt). Here, ^ab {t, x) is a symmetric rank 
2 tensor i.e., a time-dependent 2-metric such that the horizon may by viewed as a 
2-dimensional brane. In addition, we have the generators, which are the vectors 
tangent to £. When decomposing £ = dt + v^Oa w.r.t. our coordinate system, 
they appear to have a "velocity" v"^ which can also be viewed as the velocity of 
a fluid particle on the horizon. Dab is then defined as the deformation tensor of 
the horizon geometry, namely Dab = IbcD^ = \ "^^f^, where D/dt denotes 
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the Lie derivative along i = dt + v 



<^dA- It is explicitly given by 



Dab = 



{dtlAB + v'-^dclAB + Oav^^cb + dsv'^^Ac) 
\ {dtjAB + va\b + Vb\a) 



(2.63) 



where '|' denotes a covariant derivative w.r.t. the Christoffel symbols of the 2- 
geometry 'Jab- Note the contribution from the ordinary time derivative of 7^5, 
and that from the variation of the generators of velocity v"^ along the horizon. 
It is then convenient to split the deformation tensor Dab into a tracefree part 
and a trace, i.e., Dab = <^ab + \6iab, where the tracefree part cab is the 
"shear tensor" and the trace, 9 = D^ = ^^^^ df^AB + the "expansion". 
The remaining component of the Weingarten map, namely the 2-vector Qa, is 
defined as ^Ia — u-Va^ with £.n = 1. Its physical meaning can be seen from 
looking at the BH angular momentum Jh- 

Indeed, from the definition above of Jh, one finds that the total BH angu- 
lar momentum is the projection of Ha on the direction of the rotational Killing 
vector m = introduced at the begiiming of this section, so that we have 



where m^^riA is the s^s-component of Ua. It is therefore natural to define, for a 



which is similar to the result above giving the BH electric charge as the surface 
integral of the "charge surface density" aH- 

Having so defined some (fictitious) "hydrodynamical" quantities on the sur- 
face of a BH (fluid velocity, linear momentum density, shear tensor, expansion 
rate, etc.), let us now see what evolution equations they satisfy as a consequence 
of Einstein's equations. By contracting Einstein's equations with the normal to 
the horizon, we can relate the quantities just defined to the flux of the energy- 
momentum tensor T^i, into the horizon. For instance, by projecting Einstein's 
equations along £^6% one finds 




(2.64) 



BH, a "surface density of linear momentum" as tta 
With this definition, one has 



= -^^A = -^n.VAi. 




(2.65) 




(2.66) 
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where 

Dtta 
dt 



C^AB ^ {daAB + VA\B + Vb\a) ~ ^&1AB, (2.67) 



vfA 



correspond to a convective derivative, a shear and an expansion rate respectively. 
Let us recall that the Navier-Stokes equation for a viscous fluid reads 

{dt + e)TT, + v''7:,^k = -g-p + 2rjaf^k + Cd.r + h, (2.68) 

where tt; is the momentum density, p the pressure, -q the shear viscosity, cTy = 
\ + ~ Trace, the shear tensor, Q the bulk viscosity, 9 — u\ the expan- 
sion rate, and /, the external force density. The two equations are remarkably 
similar This suggests that a BH can be viewed as a (non-relativisticQl brane 
with (positiv^ surface pressure p ~ external force-density = —(-^T^a 
which corresponds to the flow of external linear momentum, surface shear vis- 
cosity rj = +y|^, and surface bulk viscosity C = ^ (in units where G = 1). 
Note, finally, that both the surface shear viscosity and the surface bulk viscosity 
apply to any type of deformed non-stationary BH. 

2.4. Irreversible thermodynamics of black holes 

In previous sections, we have introduced some electrodynamic and fluid dynami- 
cal quantities associated to a kind of dissipative dynamics of BH horizons. In ad- 
dition, following Bekenstein, we would like to endow a BH with a surface density 
of entropy equal to a dimensionless constant a (in units where h ^ G ~ c — 1). 
Any dissipative system verifying Ohm's law and the Navier-Stokes equation 
is also expected to satisfy corresponding thermodynamic dissipative equations, 
namely Joule's law and the usual expression of the viscous heat rate proportional 
to the sum of the squares of the shear tensor and of the expansion rate. More pre- 
cisely, we would expect to have a "heat production rate" in each surface element 
dA of the form 



q = AA 



(2.69) 



■^The non-relativistic character of the BH hydrodynamical-like equations may seem surprising in 
view of the "ultra-relativistic" nature of a BH. This non-relativistic-looking character is due to our 
use of an adapted "light-cone frame" (£, n,eA)- It is well-known that light-cone-gauge results have 
a distinct "non-relativistic" flavour. 

^This is consistent with the idea that the BH surface pressure must counteract the self-gravity. 
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where p is the surface resistivity, and r] and C the shear and bulk viscosities. In 

addition, one expects that this heat production rate should be associated with a 
corresponding local increase of the entropy s = adA contained in any local 
surface element of the form 

^ = ^. (2.70) 



a 

Sttq ' 



with a local temperature T expected to be equal to T 

Remarkably, the "scalar" (i'^i") projection of Einstein's equations, (i.e., , the 
Raychauduri equation) yields an evolution law for the entropy s = adA of a 
local surface element which is very analogous to what one would expect. Indeed, 
it yields 



ds d^s _ dA 



(2.71) 



where T = where 2r]aABcr + QO are exactly the expected viscous 
contributions, and where p{k, — anv^ is Joule's law. 

The only unexpected term in this result is the second term on the l.h.s.; this term 
goes beyond usual near-equilibrium thermodynamics (which involves only the 
first order time derivative of the entropy), and is proportional to the second time 
derivative of the entropy density and to a time scale r = ^ . It is interesting to note 
that, for the value a = 1/4, corresponding to the Bekenstein-Hawking entropy 
density, t is equal to i.e., the inverse of the lowest "Matsubara frequency" 
associated to the temperature T (one also notes that t ^ D = 21), where D, 
V are the diffusion constants of [14]). The minus sign in front of this new term 
is also a particularity of BH physics. In the approximation of a constant r, and 
in solving for the non-equilibrium second law of thermodynamics, one finds that 
the rate of increase of entropy is given by 

" f°''*'e-^a)in. (2.72, 



dt Jt T \T 

i.e., it is defined not as the value of the heat dissipated instantaneously, nor as 
an integral over the past heat dissipation, but as an integral over the future. This 
highlights the acausal nature of BHs, i.e., a BH is defined as a null hypersurface 
which will become stationary in the far future . As such, it has to anticipate any 
external perturbation. Failing this, the null hypersurface would generically tend 
either to collapse, or blow up toward oo. 



We also note that the ratio of the shear viscosity r] = l/(167r) to the entropy 
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density s — s/ dA = d is given by 

r] _ 1 1 
s al67r Att 



(2.73) 



where, in the second equality, we have used the Bekenstein-Hawking value for 
s = s/dA — a = J. It is interesting to note that the result f = is indeed the 
ratio found by Kovtun, Son and Starinets in the gravity duals of strongly coupled 
gauge theories [13, 14]. 



Finally, let us note another remarkable agreement between BH dissipative dy- 
namics and a rather general property of ordinary (near-equilibrium) irreversible 
thermodynamics. This agreement concerns what Prigogine has called the "mini- 
mum entropy production principle". Let us consider the total "dissipation func- 
tion" 

<j> q (2.74) 

as a functional of the velocity field {x^^x^) and of the electric potential 
(f) (a;^, x'^) in the presence of given external influences such as an external mag- 
netic field or tidal forces acting on the rotating BH. Then, -D[0] or (impos- 
ing Vyp^ = as a constraint), reach a minimum when (and only when) the lowest 

order (Einstein-Maxwell) dynamical equations for (/> or are satisfied. 



2.5. Hawking Radiation 

In this section we discuss the phenomenon of Hawking radiation, first obtained in 
Ref. [17] (we shall follow here the derivation of Ref. [18]). For simplicity, we will 
consider a 3h-1 dimensional spherically symmetric BH. We remind the reader that 
the coefficient A (r), associated to the time coordinate (here denoted as T), goes 
to zero on the horizon, so that the horizon is an infinite redshift surface. It is also a 
Killing horizon, i.e., the (suitably normalized) normal vector i = d/dt + ild/dcj) 
is a Killing vector. These points will be crucial in the following. Sinc e th e 
coefficient of the radial coordinate is defined by B (r) = -jjy^ (see Section 2.1 1, 
it is singular on the horizon. To get a good coordinate system on the horizon, we 
first factorize A{r), 

ds^ = -A(r)dT2 + ^+rW 

A{r) 

(2.75) 

= -A{r){dT'' 
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and then introduce a new radial coordinate, the so-called tortoise coordinate, 
defined by 

dr 

(2.76) 



Airy 

Note that as r ^ r-|_, A (r) ~ (^)^ ^ where (^)j, is (as mentioned 
above) equal to twice the surface gravity g. This implies 

dr 



''+)(Sr)r+ (2.77) 

In (r — r_|_) 

such that as r ^ +oo, ~ r + 2Mlnr and as r — s- r+, ~ . The 

line element can thus be re-written as 

ds^ = ~A (dT - dr,) (dT + dr,) + r^df^^. (2.78) 

We now switch to the so-called Eddington-Finkenstein coordinates, (t, r, 0, (p), 
where the combination t = T+r, of T and r^, remains regular across the (future) 
horizon 1^ and define t = T + r*. The time translation Killing vector d/dt 
coincides with the usual one d/dT. In terms of these new coordinates the metric 
reads 

ds^ = -A{r)dt^ + 2dtdr + r^dil^. (2.79) 

This metric is now regular (with a non-vanishing determinant) as the radial coor- 
dinate r penetrates within the horizon, i.e., becomes smaller than r^. 

Having defined a regular coordinate system, we now consider a massless 
scalar field, the dynamics of which is given by the massless Klein Gordon equa- 
tion 

o^Dgip^ -^d^ i^gg'^'d.^) (2.80) 

The solutions to this equation in a spherically symmetric and time-independent 
background are given by mode functions which are themselves given simply by 
products of a Fourier decomposition into frequencies, spherical harmonics and a 
radial dependence and thus read 

^^irn (T, r, 0, ^) = 4^ """"" ^""^ ^^"^ ^) ■ (2.81) 
v/27r|w| r 



^Given that the horizon is an infinite redshift surface, it takes an infinite time T to fall into it, while 
~ In (r — r+) /2g goes to — oo at the horizon. The sum of the two remains, however, finite. 
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The problem is then reduced to solving a radial equation with the radial coordi- 
nate r*, 

'^^+y-Ve{r{r,))]u = 0. (2.82) 

In the case of the Schwarzschild metric (i.e., when A{r) = 1 — 2M/r) the effec- 
tive radial potential Vt (r) is given by 

= (^1 - — j (^^^ j . (2.83) 

An essential point to note is that the effective potential vanishes both at oo (like 
a massless centrifugal potential £{£ + l)/r^), and at the horizon (where it is 
proportional to A (r)). Therefore, in these two regimes (which correspond to 
r» — > ±cxd), the solution of the wave equation behaves essentially as in flat space 
(see FIG. [3]l. The effect of the coupling to curvature is non-negligible only in 
an intermediate region, where the combined effect of curvature and centrifugal 
effects yield a positive potential barrier In turn, this potential barrier yields a grey 
body factor which diminishes the amplitude of the quantum modes considered 
below, i.e., those generated near the horizon and which must penetrate through 
the potential barrier on their way towards oo. Far from the potential barrier, the 
general solution for (p is 

(fiuj^m ~ , , , -Ylm {0, ■ (2.84) 

The quantification of the scalar field ip is rather standard (and similar to what one 
does when studying the amplification of quantum fluctuations during cosmologi- 
cal inflation). The quantum operator for the scalar field is decomposed into mode 
functions, i.e., the eigenfunctions of the Klein Gordon equation, with coefficients 
given by creation and annihilation operators. The subtlety lies, however, in the 
definition of positive and negative frequencies. 

Let us start by formally considering the simpler case of a quantum scalar field 
'f{x) in a background spacetime which becomes stationary both in the infinite 
past, and in the infinite future. In that case, one can define positive and negative 
frequencies in the usual way, in the two asymptotic regions t zLoo. The coeffi- 
cient of the positive frequency modes, say p{x), then defines an annihilation op- 
erator a. However, there are two sorts of positive-frequency {p{x)), and negative- 
frequency in{x)) modes. The in ones pf^ (defined in the asymptotic region 
t — oo, and then extended everywhere by solving the Klein-Gordon equa- 
tion), and the out ones n°"* (defined in the asymptotic region t +oo). 
The operator-valued coefficients of these modes define some corresponding in 
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and out annihilation or creation operators, so that we can write the field operator 
tf{x) as (here j is a label that runs over a basis of modes) 



E 



(2.85) 



Here, the modes are normalized as {pi,pj) = 5ij and {rii^rij) = — (5^, where 
'( , )' is the Klein-Gordon scalar product, i.e., ~ i J da^^ {ip\df^Lp2 — d^ip*i'^2)- 
The operators and are correspondingly normalized in the usual way as 

Oi, = 6ij. One defines both an in vacuum \in) and an out vacuum \out) as 

the states that are respectively annihilated by a*" or a°"*. Then, the phenomenon 
of particle creation corresponds to the fact that the out vacuum differs from the 
in one. More quantitatively, the expectation value of the number of out particles, 
in the mode labelled by i, which will be observed when the quantum field is in 
the in vacuum state is given by 

(iV,) = {in\{a'r')^ ar'\in) 

|2 (2.86) 



where we have introduced the transition amplitude Tij = , n™) from an ini- 
tial negative frequency mode n*" into a final positive frequency one These 
transition amplitudes (also called Bogoliubov coefficients) enter the calculation 
because, as is easily deduced from the double expansion of the field (p {x) above, 
they give the part of a°"* which is proportional to (a™)^- The application of the 
previous general formalism to the BH case is delicate since a BH background 
is not asymptotically stationary in the infinite future (because of the BH interior 
where the Killing vector d/dt is spacelike), and is asymptotically stationary in 
the infinite past only if we do consider explicitly the collapse leading to the for- 
mation of a BH from an initially stationary star However, Hawking showed how 
to essentially bypass these difficulties by focussing on two types of modes: 

• The high-frequency modes coming from the infinite past, which reach the 
horizon with practically no changes (because of their high-frequency nature) and, 

• the outgoing modes, viewed in the asymptotically flat region and in the far 
future. 

Concerning the outgoing modes, they can be unambiguously decomposed in 
positive- and negative-frequency parts, because, as explained above, their asymp- 
totic behaviour is given by a sum of essentially flat-spacetime modes, ( |2.84| i. One 
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HORIZON 



Figure 3. This figure is a scliematic representation of the effective gravitational potential in the 
neighborhood of a BH. Note that as far as the particles are concerned, the spacetime is essentially 
flat both at infinity and near the horizon. The tidal-centrifugal barrier that separates the horizon from 
infinity gives rise to the grey body factor. 

then defines the outgoing ^^''s as being proportional to e~'"^^~''*) with a pos- 
itive UJ. 

Let us now focus on the definition of positive- and negative-frequency modes 
near the horizon. We recall that, as mentioned at the beginning of this sec- 
tion, there is a physically infinite redshift between the surface of the horizon 
and asymptotically flat space at infinity. If one is interested in particle creation 
with a finite given frequency, as observed at infinity, the corresponding wave 
packets will have very high frequency near the horizon and can therefore be ap- 
proximated by very localized wave packets. Given that the spacetime geometry 
in the vicinity of the horizon is regular, with a finite radius of curvature, it can 
be regarded as a piece of flat spacetime locally if one looks in a small enough 
region. In this approximation, the calculation can be performed in a single step. 

We wish to compute the average number of final outgoing particleaj seen in 
the in vacuum. Then the average number of outgoing particles of type p""* is 

'Note that the "out" label, in the general discussion of particle creation above, referred to "final" 
particles (as defined in the final, asymptotic, stationary spacetime background). In the case of a 
BH background, the "final" spacetime is made of two separate asymptotic regions: (i) the outgoing 
wave region at spatial infinity, and (ii) the vicinity of the (spacelike?) singularity within the BH. 
The definition of positive- and negative-frequency modes in the latter region is ill-defined. However, 
luckily, the calculation of the physically relevant flux of final, outgoing modes can be performed 
without worrying about the physics near the BH interior singularity. In other words, it is enough to 
consider as "out" positive-frequency modes p°"* only the ones outgoing at spatial infinity (i.e., on 
"scri+"), though they do not constitute a complete basis of final modes. 
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given by \Tij\^, where Tij = [p°'"*,nj^) is the transition amplitude from 
an initial negative frequency mode into a final outgoing positive frequency 
one p^"* (recorded at spatial infinity). To compute this transition amphtude, we 
need to describe what is an initial negative frequency mode n^". As said above, 
Hawking suggested that only high-frequency initial modes are important, and that 
they essentially look the same (some kind of WKB wave) in the real in region (in 
the far past, before the formation of the BH) as in the vicinity of the horizon. Our 
technical problem is then reduced to characterizing what is a negative frequency 
mode as seen in a small neighborhood of the horizon, which looks hke the 
Minkowski vacuum. 

To do this, it is convenient to have a technical criterion for characterizing 
positive and negative frequency modes in (a local) Minkowski spacetime. Lo- 
cally, one can perform a Fourier decomposition of the wave packet and use 
the mathematical fact that the Fourier space properties are mapped onto ana- 
lytic continuation properties in x-space. This relation can then be used to define 
positive and negative frequency modes. This is easy to see. Consider a gen- 
eral negative frequency wave packet in Minkowski spacetime. It has the form 
(p- (x) ~ d^k(p (k) e^''"^" where fc^ is timelike -or-null and past-directed, 
i.e., E C^. We now perform a complex shift of the spacetime coordi- 
nate, Xij, ^ x'^ + iy'^, where is timelike-or-null and future-directed (i.e., 
yij- g ^7+)^ (jien^ the e**'"^" term will be suppressed by a e"*^""" term, where the 
scalar product k^y^ is positive because it involves two timelike vectors that point 
in opposite directions (we use the "mostly plus" signature). This ensures that a 
negative-frequency wavepacket can indeed be analytically continued to complex 
spacetime points of the form -\- iy^, with y^ e C+. 

The strategy for applying this criterion to characterizing negative-frequency 
modes nj^ in the vicinity of the horizon is then the following. One starts from 
a wavepacket which is not purely a "negative-frequency" one near the hori- 
zon, but which has the property of evolving into an outgoing positive-frequency 
wave packet (so that it will have a non-zero transition amplitude to some p°"*). 
Then, one modifies the initial wavepacket so that it becomes a purely negative- 
frequency mode n" near the horizon. 

When looking at a wavepacket of the form of 95°"* {t, r) oc g-^'^C^-''*) (with 
positive uj) just outside the horizon, one must first switch to well-defined coordi- 
nates to examine its physical content. We therefore replace the Schwarzschild- 
type time coordinate T by the Eddington-Fikenstein time coordinate, t = T -|- r* 
which is regular on the horizon. After rearranging terms according to T — r* = 
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(T + r*) — 2r* the previous (outgoing, positive-frequency) wave reads 

K"* r)] oc e-^'^C^-'^*' = e-^'^t^2iu,r, 

= e"*'^*e'9 '"^''"''+^ (2.87) 
= e-''^*(r-r+)^ . 

This describes the behaviour, just outside the horizon, of a wavepacket which 
will become (modulo some grey-body factor) an outgoing positive frequency 
wavepacket at oo. However, locally on the horizon, it is neither a positive nor 
a negative frequency wavepacket because, at this stage, it is defined only out- 
side the horizon, but not inside. Let us now show how one must continue this 
wavepacket inside the horizon, so that it becomes a genuine negative-frequency 
wavepacket "straddling" the horizon. Using the criterion explained above, we can 
"continue" the wavepacket inside the horizon by a suitable analytic continuation. 
More precisely, we need an analytic continuation of the form x^^ + iy^, 

where belongs to the future lightcone to ensure that we shall then be dealing 
with a local negative frequency wavepacket. It is easy to see, from a spacetime 
diagram of the Ughtcone on the horizon, that the vector d/dr is everywhere nuU 
and past-directed, such that r ^ r — e, where e > 0, is everywhere null and fu- 
ture directed. The analytic continuation of {t,r) tor ^ r — ie will therefore 
define for us a good local negative-frequency mode n*" = n^^. One easily sees 
that this analytic continuation in r generates a new component to the wavepacket 
which is located inside the BH (i.e., for r < r+). More precisely, a one line 
calculation yields 

n^em{r,t) = N^^^f {t,r - is) 



= N^[0ir~r+)<p°J'{r-r+)+ (2.88) 
eT^(r+-r)^r (r+ - r) 

where the second term is the wavefunction inside the horizon that has acquired 
an additional exponential factor due to the rotation e^™ in the complex plane 
from r > r+ to r < r+. The overall factor Af„ is a normalization factor (needed 
because we have extended the mode inside the BH), such that 

(r, t) T^:,,,^, (r, t)) = 8{uj- uj') ht' ^mm- , (2.89) 

from which we obtain (when remembering that was correctly normalized) 
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The physical meaning of equation ( |2.88| l is the description of the splitting of the 
in mode nj" = into a positive-frequency wave outgoing from the hori- 

zon and a wave falling from the horizon towards the singularity. One can read 
off from it the needed transition amplitude Tij ~ (p°"*,nj"). It is essentially 
given by the factor N^^, which must, however, be corrected by a grey-body factor 
y^Ti (uj) taking into account the attenuation of the outgoing wave e^*'^'^-'""''*^ as 
it crosses the curvature + centrifugal potential barrier Ve (r) on its way from the 
horizon to oo. Then (using Fermi's golden rule), one easily finds that the general 
result ( |2.86| l yields a rate of particle creation given by 



dt ^ / 27r 



e 9 — 1 



(2.91) 



One recognizes here a thermal (Planck) spectrum (corrected by a grey-body fac- 
tor). From the Planck factor, one reads off the Hawking temperature, T — 
This result fixes the dimensionless coefficient a in the Bekenstein entropy to the 
famous result a—\, i.e., 

A 



Let us end by two final comments. First, the generalisation of the Hawking 
radiation to a more general rotating and/or charged BH is given essentially by 
replacing in the result above the frequency w by w — wq where ujq — m£L + eV 
exhibits the couplings of the created particles to the angular velocity and the 
electric potential V of the BH. Then, in astrophysically realistic conditions, the 
"Hawking" part of the particle creation (i.e., the thermal aspect) is too small to 
be relevant, while the combined effect of the grey body factor and of the zero- 

temperature limit of (e s —1) yield potentially relevant particle creation 
phenomena in Kerr-Newman BHs, associated to the "superradiance" of modes 
with frequencies [l < uj < cjq, where /i is the mass of the created particle (see, 
e.g., [18] for more details and references). We conclude by noting that the situ- 
ation just described is not only technically similar to the one in the inflationary 
scenario for cosmological perturbations, but also physically similar in that, in 
both cases, transplanckian frequency modes in the ultraviolet are redshifted to a 
finite, observable frequency. 



3. Experimental tests of gravity 

Before discussing various possibilities of string-inspired phenomenology (and of 
possible string-inspired deviations from Einstein's theory of General Relativity) 
we give an overview of what is known experimentally about the gravitational 
sector 
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3.1. Universal coupling of matter to gravity 

The standard model of gravity is Einstein's General Relativity (GR). In GR, aU 

fields of the standard model of particle physics (SM) are universally coupled to 
gravity by replacing the flat spacetime metric 1]^^ by a cmved spacetime one . 
In "standard GR" one also assumes that gravity is the only long range coupling 
(apart from electromagnetism). We shall see below, how the presence of other 
long range interactions (coupled to bulk matter) modify the usual "pure GR" 
phenomenology. The action for the matter sector, Ssm, has the structure 

^-^g^^'^D^HD^H - ^gV {H) 

(3.1) 

where D denotes a (gauge and gravity) covariant derivative, while the dynamics 
of gfj^i, is described by the Einstein-Hilbert action, Seh, 

/4 
d^x^^^g'^'R^^ig). (3.2) 

The total action is therefore given by 

S = Seh [9^,.] + SsM A^, H; g^,] , (3.3) 

and its variation w.r.t. yields the well-known Einstein field equations 

1 SttG 
Rfiv — 2^9iiv = —^T^j, (3.4) 

where T'"^ = ^fg^^ ■ The universal couphng of any type of particle to 
5^1, is made manifest in Ssm while Seh contains all the information on the 
propagation of gravity. For instance, expanding Seh in powers of (where 
g^v = Vnu + /^jui/)) one obtains, at quadratic order in /i^^,, the spin-two Pauli- 
Fierz Lagrangian. Higher orders in /i^,^ contain an infinite series of nonlinear 
self-couplings of gravity: ddhhh, ddhhhh, etc. As we shall see, this nonlinear 
structure has been verified experimentally to high accuracy (both in the weak- 
field regime, where the cubic vertex ddhhh has been checked, and in the strong- 
field regime of binary pulsars, where the fully nonlinear GR dynamics has been 
confirmed). In the following, we discuss, successively, (i) the experimental tests 
of the couphng of matter to gravity, and (ii) the tests of the dynamics of the grav- 
itational field: kinetic terms (describing the propagation of gravity), and cubic 
and higher gravitational vertices. 
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The universal nature of matter's coupling to gravity, i.e., the coupling of matter 
to a universal deformation of spacetime, has many experimental consequences. 
These experimental consequences can be derived by using a simple theorem by 
Fermi and Cartan. Given any pseudo-Riemannian manifold, for instance a curved 
spacetime endowed with a metric g^j^, and given any worldline C in this space- 
time (not assumed to be a geodesic), there always exists a coordinate system such 
that, all along C, g^^ [x^^ = r/^j^ + O (x^), where x denotes the spatial devi- 
ation away from the central worldline C It is important to note that there is no 
linear term in x, but only x'^ effects, i.e., tidal effects. There exists a very simple 
and intuitive demonstration of this Fermi-Cartan theorem. Let us view the curved 
manifold as being some "brane" embedded within a flat ambient auxiliary man- 
ifold. For instance, consider an ordinary 2-surface E within a three-dimensional 
flat euclidean space. Given any (smooth) curve C traced on S, we can take a 
flat sheet of paper and progressively "apply" (or "fit") this sheet on E along the 
curve C. The orthogonal projection of S onto this applied flat sheet defines a 
map from S to a coordinatized flat manifold which has the property enunciated 
above. Note that, in this "development" of the neighbourhood of C within E onto 
a flat sheet, the shape (as seen on the flat sheet) of the "developped" curve C is 
generically not a straight line. It is only when C was a geodesic line on S, that its 
development will be a straight line. This proof, and its consequences, are valid 
in any dimension and signature. 

Here, we have in mind applying this result to the "center of mass" worldline C 
of an arbitrary body moving in a background spacetime, or more generally of any 
sufficiently small laboratory (containing several bodies, between which we can 
neglect gravitational effects). We assume that we can neglect the backreaction of 
the body (or bodies) on the spacetime. In the approximation where we can neglect 
the tidal effects (linked to the O (x^) terms in in Fermi coordinates), we can 
consider that we have a body, or a smaU lab, moving in a flat spacetime. In other 
words, the theorem of Fermi and Cartan tells us that we can essentially "efface" 
(modulo small, controllable tidal effects) the background gravitational field 17^1, 
all along the history of a small lab, or a body. This "effacement property" is 
telling us, for instance, that the physical properties we can measure in a small lab 
will be independent of where the lab is, and when the measurements are made. In 
particular, all the (dimensionless) coupling constant^that enter the interpretation 
of local experiments (such as various mass ratios, the fine-structure constant, etc.) 
must be independent of where and when they are (locally) measured (constancy 
of the constants). A second consequence of this effacement property is that local 
physics should be Lorentz 5*0(3, 1) invariant, because this is a symmetry of 



**We assume here that the cutoif length scale t = 1/A of any low-energy effective QFT descrip- 
tion of the physics in a small lab is is fixed, when measured in units of ds = ^g^^daj'^dx". 
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the (approximate) flat spacetime appearing after one has effaced the tidal effects 
{local Lorentz invariance). 

Moreover, in absence of coupling to other long-range fields (such as electro- 
magnetism for a charged body), the center of mass of an isolated body (viewed 
as moving in a flat spacetime) must follow a straight worldline (principle of iner- 
tia). We therefore conclude (by the theorem above) that C has to be a geodesic in 
the original curved spacetime. This is true independently of the intemal proper- 
ties of the object. One may thus conclude that isolated neutral bodies fall along 
geodesies independently of the internal properties of the object, since at no point 
in the demonstration had we to rely on any internal properties of the object. This 
is therefore a proof of the weak equivalence principle, i.e., all bodies in a gravi- 
tational field fall with the same acceleration. Note, once again, that the absence 
of other long range fields besides gij,^ that could influence the object considered 
is crucial. 

Finally, another universality property, that of the gravitational redshift, may 
be shown by a comparison of the GR formulation with the Newtonian one. In 
lowest order approximation, the deviation of the goo component of the metric 
from ?7oo = — 1 is twice the Newtonian potential U{x). Indeed, comparing the 
action for a geodesic, 

SE = -m j dt^-g^^x^'x'' (3.5) 

with 



'Newton 



-I 



^ mx^ + mil {x) 



one finds goo = -1 + {x) + O (^) where U = 



E 



(3.6) 

Grua 



Experimentally, one may transfer electromagnetic signals from one clock to 
another identical clock located in a gravitational field. If we are in a stationary 
situation (i.e., if there exists a coordinate system w.r.t. which the physics is in- 
dependent of time x^ = ct), the time translation invariance of the background 
shows that electromagnetic signals will take a constant coordinate time to prop- 
agate from clock 1 to clock 2. We can then use the hnk drj = \/—goo{xi)dti 
between the proper time (at the location of clock t (z = 1, 2)) and the correspond- 
ing coordinate time, as well as the (approximate) result above for goo- Finally, 
we conclude that two identically constructed clocks located at two different po- 
sitions in a static external Newtonian potential exhibit, when intercompared by 
electromagnetic signals, the (apparent) difference in clock rate 
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This gravitational redshift effect is proportional to the difference in the Newto- 
nian potential between the two locations, independently of the constitution of the 
clocks (say Hydrogen maser, or Cesium clock, etc.). This is a property known as 
the universality of the gravitational redshift. 

The various consequences, discussed above, of the universal character of the 
coupling of matter to gravity are usually summarized under the generic name of 
equivalence principle. In the next section, we discuss the experimental tests of 
the equivalence principle and their accuracy. 



3.2. Experimental tests of the coupling of matter to gravity 
3.2.1. How constant are the constants? 

The best tests of the "constancy of the constants" concern the fine structure con- 
stant a = /he ~ 1/137.037, and the ratio of the electron mass to that of the 
proton ^ (see Ref. [21] for a review). There exist several types of tests, based, 
for instance, on geological data {e.g., measurements made on the nuclear decay 
products of old meteorites), or on measurements (of astronomical origin) of the 
fine structure of absorption and emission spectra of distant atoms, as, e.g., the 
absorption lines of atoms on the line-of-sight of quasars at high redshift. Such 
kinds of tests all depend on the value of a. There exist, in addition, several lab- 
oratory tests such as, for example, comparisons made between several different 
high- stability clocks. However, the best measurement of the constancy of a to 
date is the Oklo phenomenorj^ It sets the following (conservative) limits on the 
variation of a over a period of two billion years [22-24] 

Q,Oklo _ today 

-0.9xl0-^< < 1.2x10"^. (3.8) 

Q^today ^ ^ 

Converting this result into an average time variation, one finds 

- 6.7 X 10-^^yr"i < - < 5 x IQ-^Vr"^- (3-9) 

a 

' The Oklo phenomenon was discovered by scientists at the Commissariat a I'energie atomique 
(CEA) in France. A study of the uranium ore in a Gabonese mine revealed an unusual depletion in 
jj2'i5 (yggfj flssion reactors) w.r.t. the usual proportion. Uranium ore is a mix of two isotopes, with, 
in usual samples, 99.28% (7^^** and 0.72% C/^^^. By contrast, the Oklo ore had only < 0.72% of 
was realized that a natural fission process took place, prompted by the presence of ground 
water, in Oklo some two billion year's ago, and lasted for about two million years. Scientists analysed 
in detail the 2 billion year-old fission decay products. One can then infer from these measurements 
the scattering cross-sections of slow neutrons on various isotropes. Then, modulo some further as- 
sumptions about the dependence of various nuclear quantities on a, one could constrain the variation 
of a between the time of the fission reaction (roughly two billion years ago) and now. For details 
about the analysis and interpretation of Oklo data see [22] and references therein. 
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Note that this variation is a factor of ^ 10^ smaller than the Hubble scale, which 
is itself ~ 10^^° yr^^. Comparably stringent limits were obtained using the 
Rhenium 187 to Osmium 187 ratio in meteorites [25] yielding an upper bound 
^ = (8 ± 8) X 10~^ over 4.6 x 10^ years. Laboratory limits were also obtained 
from the comparison, over time, of stable atomic clocks. More precisely, given 
that - a for electrons in the first Bohr orbit, direct measurements of the vari- 

c 

ation of a over time can be made by comparing the frequencies of atomic clocks 
that rely on different atomic transitions. The upper bound on the variation of a 
using such methods is ^ = (—0.9 ± 2.9) x 10^^^ yr^^ [26]. It should be men- 
tioned that a few years ago claims were made concerning observational evidence 
of non-zero time variations of a and ^ from analyses of some astronomical 
spectra (see Ref. [21]). Other recent astronomical data indicate no variability 
of these constants (see Ref. [21] and the chapter 18 of the Review of Particle 



3.2.2. Tests of local Lorentz invariance 

We should first mention that the Michelson-Morley experimenj^ has been re- 
peated (with high accuracy) and strong limits have been obtained on a possible 
anisotropy of the propagation of light. In its modem realizations (Brillet and 
Hall, 1979), it has been performed with laser technology on rotating platforms. 
This experiment is now viewed as a test of the isotropy of space on the moving 
Earth, and thereby as a test of local Lorentz invariance. There also exists another 
idea for testing the isotropy of space, and although its interpretation is not totally 
clear, it is a conceptually interesting idea. This is why we choose to outline it in 
these lectures. 

For simplicity, consider the hydrogen atom. Assuming the isotropy of space, 
i.e., the existence of a SO(3) symmetry, we know that there should exist a degen- 
eracy in the energy levels, given by the magnetic quantum number m. However, 
it is interesting to understand how the SO(3) symmetry comes about dynamically 
(and therefore, how it might be dynamically violated). The Hamiltonian for the 
electron is given by 



where the first term is the kinetic term (A being the Laplacian), and the second 
term is the Coulomb potential. Note that in fact, A = and — SijX^x^ , 

such that both terms depend on the same spatial structure 5ij, the flat metric, 

'"Available on http://pdg.lbl.gov/ 

"First performed as part of a series of experiments, beginning in Potsdam in 1881 (by Michelson 
alone) and then in the US until 1887 (by both Michelson and Morley) to test the existence of the 
aether. 




for references). 
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thereby ensuring the SO(3) symmetry. However, both terms also come from 
an underlying field theoretic formulation: (i) the non-relativistic electron kinetic 
energy term oc A = S^^dij comes from the kinetic term in the Dirac action, 
'ip'yi^dfiip — rmpip, with {"f^,^"} = rj^i,, while (ii) the e^/r term is the static 
Green's function of the electromagnetic field, which comes from inverting the 
kinetic term of the photon ri°'^ri^'' FapF^i,, which manifestly depends, by as- 
sumption, on the same spacetime metric 77^^. Einstein assumed that, in order 
to take into account the coupling to gravity, it was sufficient to replace r;^^ by 
the same g^jy both for the electron and the photon. By contrast, let us consider 
the possibility that electrons ("matter") and photons ("electromagnetism") have 
a different coupling to gravity, e.g., described by saying that they couple to two 
different (spatial) metrics, say 



Then, computing the new propagators for the electron and the photon in their re- 
spective metrics, one finds that the SO(3) symmetry would be violated by tensor 
terms, appearing in the Hamiltonian, of the form 5H ~ ^hij^^. This is a 
violation, at a deep level, of the universality discussed in the previous section. 
The usual SO(3) symmetry imphes that all energy levels with magnetic quantum 
number m are degenerate. But if tensor terms violating SO(3) were to exist, then, 
observables effects would include potentially measurable quadrupole-type split- 
tings in the energy levels, which, apphed to the atomic nucleus (whose energy 
levels are a more sensitive probe of anisotropy), are oc {I M\Qij\I M), where 
/ and M are the nuclear spin quantum numbers, and where Qij is a symmetric 
tracefree tensor operator that couples to the tracefree part of hij . Such types of 
measurements have been performed on the energy levels of nuclei with impres- 
sively high accuracy, the current upper bound being 



The universahty of space is thus valid to one part in 10 ', showing how delicate 
Einstein's postulate is. 

3.2.3. Universality of free fall 

The most recent limits on the deviation from the universality of free fall have 
been obtained by Eric's Adelberger's group [27]. In particular, they compared the 
acceleration of a Beryllium mass and a Copper one in the Earth's gravitational 
field and found 




Sii < 10 



-27 



(3.12) 




Be-Cu 



= (-1.9 ±2.5) X 10 



,-12 



(3.13) 
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where Aa = — o-Cu- Other limits exist, such as, for instance, the fractional 
difference in acceleration of earth-core-hke (~ iron) and moon-mantle-Uke (sil- 
ica) bodies, 



Aa 

a 



= (3.6 ± 5.0) X 10"^^ (3.14) 

Earth— core— Moon— mantle 



There are also excellent limits concerning celestial bodies. In particular the 
possible difference in the accelerations of the Earth and the Moon towards the 
Sun have been measured using laser ranging (with 5 mm accuracy) with retro- 
reflectors (corner cubes) placed on the Moon, giving the result [28] 



Aa 



j = (-1.0 ±1.4) X 10"^^. (3.15) 

/ Earth-Moon 



One should, however, remember that only a fraction ( ~ 1/3) of the Earth mass 
is made of iron, while the rest is mostly sihca (which is the main material the 
Moon is made of). As, independently of the equivalence principle, silica must 
fall like silica, one looses a factor 3, so that the resulting bound on a possible 
violation of the equivalence principle is only around the 5 x 10 level, which 
is comparable to laboratory boimds. 

3.2.4. Universality of the gravitational redshift 

We conclude the section on the tests of the coupling of matter to gravity by just 
mentioning that the universality of the gravitational redshift, namely the apparent 
change in the frequencies of two similar clocks in a gravitational field, has been 
tested by comparing the frequencies of hydrogen masers at the Earth surface and 
in a rocket. Vessot and Levine (1979) in Ref. [29] verified that the fractional 
change in the measured frequencies is consistent with GR to the 10~^ level: 



^ = (i±io-)^^ 



The universality of this redshift has also been verified by measurements involving 
other types of clocks. 

3.3. Tests of the dynamics of the gravitational field 
3.3.1. Brief review of the theoretical background 

Until now we have only considered the couphng between matter and gravity, and 
various tests of its universahty. We now discuss the tests of the dynamics of 
the gravitational field, i.e., tests probing either the propagator of the gravitational 
field, or the cubic or higher order gravitational vertices (for more detailed reviews 
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see Refs. [30, 31]). We first consider the weak field regime, regime in which we 
can write g^,^ = rj^i, + /i^^, where h^,^ is numerically much smaller than one. 
For instance, in the solar system, h^v ~ 10~^ on the surface of the Sun, ~ 10~^ 
on the Earth orbit around the Sun, or ^ 10~^ on the Earth surface. With values 
so small, it is clear that the solar system will not allow one to test many nonlinear 
terms in the perturbative expansion of g^^,. 

We start by considering the gravitational interaction between two particles 
of masses uia and tob. At linear order in /i^,^, we will have an interaction 
corresponding to the following (classical Feynman-like) graph 




To compute explicitly what the preceding graph means, we must start from the 
full action describing two gravitationally interacting bodies A and B: 

S = -ruA J dsA J dsB + J J^^- (3-17) 

Expanding S in the deviations of g^i, away from rj^i,, one obtains (denoting 
h = v'^'K.) 



S = —rriA 
1 
2 



+ 



0{h'^T) + 0{h^), 

(3.18) 

where T'^ is the (flat-space hmit) of the stress-energy tensor of particle A (given 
by a (^-function localized on the worldline of A), and where the kinetic term of 
/lyai/ is the one corresponding to the harmonic gauge {i.e., dy {y/gg^'^) = 0). 
Inverting this kinetic term yields for the following lowest-order equation 
(corresponding to Einstein's equations at hnearized order) 

□V = (t^m. - • (3.19) 

with T^ij = rinariufj{T2^ +T'g^). We can then "integrate out" h^^ by solving the 
latter field equation for h, and replacing the result in the original action. Modulo 
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self-interaction terms oc T^'^O ^P^upaTj^ , the action then splits into the sum 
of three terms, a term —rriA J \/—ilfj,i''^^A'^^% describing the free propagation 
of body A, a similar term for B, and an interaction term. 



■5'"' = / ^aO-' { T^. j^T^V,. ] • (3.20) 



More explicitly, if we introduce the scalar Green's function G{x), such that 
DG{x) = —A'k6^{x), this lowest-order interaction term reads 

G {XA {sa) - XB {sB))mBUBpUBa, 

in which one easily identifies the usual structure of a Feynman graph (namely 
the one depicted above), with the coupling constant G in front, and a graviton 
propagator (comprising the scalar Green's function, together with the spin-two 
projection operator P^^, which can be read off the previous explicit result) sand- 
wiched between two source terms. 

In the stationary approximation, the scalar Green's function reduces to the 
usual Newtonian propagator 1/r. If one further neglects the relative velocity 
of the two worldlines one can replace the spacetime velocities it^ and Ug by 
(1, 0, 0, 0, • • • ). This yields the usual Newtonian interaction term G J dtmArriB/r. 
However, the "one-graviton exchange" diagram above contains many Einsteinian 
effects that go beyond the Newtonian approximation. To compute them explic- 
itly, we first need the explicit expression of the relativistic scalar Green's function 
G{x). As we are deriving here the part of the gravitational interaction which is 
"conservative" (i.e., energy conserving), we must use the time-symmetric (half- 
advanced half-retarded) Green's function. In four dimensions, it is given by 
6 [{xa — xb)^)- It is the sum of two terms (a retarded and an advanced one), as 
depicted in FIG.|4] Note in passing that this classical time-symmetric propagator 
corresponds to the real part of the Feynman propagator. Indeed, for a massless 
scalar particle in 13 = 4 the Feynman propagator (in x space) is proportional to 

, ' . = iPP^ + ttS (x^) , (3.22) 
x"^ + le x^ 

where the first term on the r.h.s. is a distributional "principal part" (it is pure 
imaginary and "quantum"), while the second (real) term is the classical contribu- 
tion (classically the interaction propagates along the light cone, see FIG.|4]). Note 
that, contrary to the Newtonian picture where the interaction is instantaneous, we 
have here an interaction which depends both on the future and on the pasj*^ 
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Figure 4. Time-symmetric half-advanced half-retarded contributions to the gravitational interaction 
between particles A and B. 

When considering the case (of most importance in applications) where A and 

B move slowly relative to c, the time-symmetric propagator can be expanded 
in powers of 1/c. The first term in this expansion yields the usual Newtonian 
instantaneous interaction, while all the higher-order terms can be expressed in 
terms of successive derivatives of the positions of A and B, (so that the acausal- 
ity formally disappears, and is replaced by a dependence of the Lagrangian on 
derivatives higher than the velocities, i.e., accelerations, derivatives of accelera- 
tions, and so forth. Actually, such higher-derivative terms start appearing only at 
the so-called "second post-Newtonian" (2PN) order, i.e., the order C (^). Such 
higher-order post-Newtonian (PN) contributions are important for some appli- 
cations (binary pulsars, coalescing black holes), and have been computed up to 
the 3PN {O (^)), as well as 3.5 PN (O (^)) levels. Here we shall consider 
only the first post-Newtonian, IPN, level, i.e., O (^)- At this level, the action 
can be written entirely in terms of the velocities of A and B (taken at the same 
instant t in some Lorentz frame). By expanding the time -symmetric acausal one- 
gravi ton-exchange action written above to order one finds the following 
exphcit IPN Lagrangian (now considered for an A^-body system made of masses 



This "acausal" behaviour is due to our considering the conservative ("Fokker") action. If we were 
computing the "real" classical equations of motion of the two particles, we would use only a retarded 
Green's function. The equations of motion so obtained would then be "causal" and would automati- 
cally contain some (physically needed), time-asymmetric "radiation reaction" terms. The trick, used 
here, to employ an acausal time-symmetric Green's function is a technical shortcut allowing one to 
derive the action yielding the conservative part of the equations of motion. 
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labelled A, B = 1,2,- •• ,N.) 

^2— body 



1 Gtuatub 



2^ 



- i+ ^ 

AB L 2 
7 ^ ^ 1 _ _ ^ . „ / 1 



2c2 



(3.23) 



Note that the coefficients 3/2, 7/2, 1/2, etc., arise from the spin 2 nature of the 
graviton, i.e., they are uniquely fixed by Einstein's propagator. 

When considering a gravitationally bound AT-body system, we must remem- 
ber that there is a link ^ due to the "virial theorem". This link says 

2 

that the ^ contributions in the one-graviton-exchange graph considered above 
must be completed by computing non-Unear interaction graphs containing more 
gravitons, namely the ones of order G^/c^ involving two powers of the coupling 
constant G. These contributions correspond to the terms 0{h'^T) or 0{h^) in the 
/i-expansion of the exact Einstein action. In terms of Feyrmiam-like diagrams, 
this means the following graphs: 






Note, in particular, that these terms involve the graviton cubic vertex (whose 
structure will therefore be probed by solar-system experiments). Note also that 
some of these terms involve only two bodies (being proportional, say, to m^m^), 
while others can involve three distinct bodies (oc mArriBmc)- The full G"^ /c^ 
result (containing both two-body and three-body terms) is found to be equal to 



^3— body 



- y 

2 ^ 



G'^iriATnB'mc 
tabTacc^ 



(3.24) 



where the factor of 1/2 is a prediction from Einstein's theory. Note that the 
summation is restricted by -B ^ ^ ^ C which allows for the two-body terms 

where B = C. 

When looking at the nonlinear diagrams above one sees some "loops" made 
by graviton propagator lines closing up on a matter worldUne. This may seem 
paradoxical because we are considering here classical gravitational effects, and 
classical theory is usually thought of as involving only tree diagrams. Indeed, if 
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we replace all our "source worldlines" (drawn above as continuous worldlines) 
by separate external sources (i.e., by replacing the line describing T'^ {x) by 
separate "blobs" on which graviton propagators start or end), we see that all the 
diagrams above open up and become tree diagrams. However, the presence of 
loops in the diagrams used here do correspond to essentially some of the same 
physical effects that "quantum loops" describe. This is particularly clear for the 
diagrams below (which are included in the classical calculations) 



It is clear that these diagrams describe the self-gravity effects of a mass uia on 
itself As such, they do describe the classical hmit of quantum loops such as the 
simplest one-loop diagram 



which describes the back action of the emission and reabsorption of a graviton 
on a quantum particle. Another similarity beween "classical loops" and quantum 
ones, is that, in practice, multi-loops are associated to the presence of multiple 
integrals which are increasingly difficult to compute[^ In addition, the loop di- 
agrams depicted in the penultimate graph lead (like quantum loops) to formally 
divergent integrals. The origin of these divergences is that we have been de- 
scribing the gravitationally interacting bodies as "pointlike", i.e., mathematically 
described by a (5-function (on a worldline). There are several ways of dealing 
with this technical problem. One can complete the formal perturbative calcu- 
lations done with point-like bodies by another approximation scheme in which 
each body is locally viewed (in its own rest frame) as a weakly perturbed iso- 
lated body. The development of such a dual perturbation method [32, 33] shows 
that the non-point-like, internal structure of non-rotating compact bodies (neu- 
tron stars or black holes) will enter their translational dynamics only at the 5PN 






Note that classical diagrams must all be computed in x-space. An increasing number of loops 
signals the presence of intermediate veilices in x-space on which one must integrate. 
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level G^/c^°), corresponding to 5 loops ! Knowing this, one expects that 
the use of a gauge-invariant regularization method for treating gravitationally- 
interacting point masses should give a physically unique answer up to 5 loops 
(excluded). Using dimensional regularization, one finds that all self-gravity ef- 
fects are unambiguous and finite at IPN and 2PN [33]. Recent work has pushed 
the calculation to the 3PN (i.e., 3 loop) level. Again, one finds, either (when us- 
ing a convenient gauge) a finite answer [34], or, when using the harmonic gauge, 
an equivalent answer after renormalizing the position of the worldline used in the 
(5-function source [35]. 

3.3.2. Experimental tests in the solar system 

The IPN-level results described above are accurate enough for describing the 
gravitational dynamics in the solar system. Testing the validity of GR's descrip- 
tion of the gravitational field's dynamics is then achieved by verifying the agree- 
ment of the coefficients introduced above with experimental measurements. In 
this section, we wiU see that these GR-predicted coefficients agree with their 
experimentally measured value to better than the 10 ~^ level for the coefficients 
entering the one-graviton-exchange term, and to about the 10~^ level for the ad- 
ditional IPN multi-graviton term. There are many observables that can be used 
to test relativistic gravity in the solar system. One may use the advance of the 
perihelion of planets, the deflection, by the local curvature, of Ught reaching the 
Earth from distant stars, the additional time delays suffered by electromagnetic 
signals compared to their flat spacetime counterparts, or also, general relativis- 
tic corrections to the Moon's motion using the laser ranging technique already 
mentioned in previous sections. 

When testing Einstein's predictions it is convenient to embed OR within a 
class of alternative gravity theories. For instance, one could consider not only the 
interaction of matter with the usual Einstein (pure spin-2) graviton but also an in- 
teraction with a long-range scalar field ^p, i.e., a spin-zero massless field, coupled 
to the trace of T**"^ with strength \/Ga{ip). This leads to an additional attractive 
force, so that the effective gravitational constant measured in a Cavendish exper- 
iment is Geff = G{\ + a^). This also modifies the v"^ /c^ terms in the two-body 
action introduced above by terms proportional to a^. These modifications are of- 
ten summarized by writing the IPN-level metric generated by an A^-body system 
in the form 



where U = GeS J2a ''^a/i'a is the (effective) Newtonian potential and where 
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the two dimensionless coefficients ^ and 7 encode the two possible (Lorentz- 
invariant) deviations from GR which enter the IPN level. For an additional 
coupling to a scalar with (^-dependent coupling strength a{ip), one finds that 
the "post-Einstein" parameter 7 is given by 7 = — 2a^/(l + a^). As for the 
other "post-Einstein" parameter /? it measures a possible modification of the 
(cubic-vertex related) three-body action £3-body yyi-j^gn above, and it is given 
by ^ = + ^^0^ / {l + a^)"^ where /3 denotes the derivative of the scalar coupling 
a w.r.t. the field Lp. 

The most accurate test of GR in the solar system is the one made using the 
Cassini spacecraft by the authors of Ref. [36]. This test is (essentially) only 
sensitive to the post-Einstein parameter 7 (i.e., it depends only on the graviton 
propagator and the coupling to matter but not on nonlinear terms). This experi- 
ment used electromagnetic signals sent from the Earth to the Cassini spacecraft 
and transponded back to Earth, and monitored the ratio between the electromag- 
netic frequency v + Aj/ recorded back on Earth to the initial frequency y. This 
ratio was used to probe the change in the geometry of spacetime in the vicinity 
of the Sun as the fine of sight moved, especially when it was nearly grazing the 
Sun. The theoretical prediction for the experimental quantity measured in this 
experiment is 



where b is the impact parameter, i.e., the distance of closest approach of the sig- 
nal's trajectory to the center of the Sun. The experimental data gave the following 
result for the parameter 7 



This confirms GR (namely 7*^^ = 0) to the 0{l) x 10"^ level. 

The three-graviton vertex can be probed by considering a body, having a non- 
negligible gravitational self-binding energy, in an external gravitational field. In- 
deed, as emphasized by Nordtvedt [37, 38], the free fall acceleration of a self- 
gravitating body is, in most gravity theories (except GR) sensitive to its gravi- 
tational binding energy. For instance, the Earth and the Moon will have, in a 
general theory, a slightly different acceleration of free fall towards the Sun. The 
effect is proportional to the combination of post-Einstein parameters 4/? — 7. Lu- 
nar laser ranging data have allowed one to put a stringent upper limit on such a 
possibility [28], namely 




(3.25) 



7= (2.1 ±2.3) X 10 



1-5 



(3.26) 



4/3-7= (4.4 ±4.5) x 10 



1-4 



(3.27) 



Thus, to date, predictions of Einstein's theory in the linear (one-graviton-exchange) 
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approximation have been verified to the 10^^ level, while some of the cubically 
nonlinear aspects have been verified to the 10^'^ level. 

The tests discussed up to now concern the quasi-stationary, weak-field regime, 
as it can be probed in the solar system. We shall now discuss the tests obtained 
in binary pulsar data, which have gone beyond the solar-system tests in probing 
part of the strongly nonlinear regime of gravity. 

3.3.3. Objects with strong self-gravity: binary pulsars 

Binary pulsars were discovered by Hulse and Taylor in 1974. Such systems are 
made of two objects going around each other in very elliptical orbits. Both ob- 
jects are neutron star^^ of which one is a pulsar, i.e., a rotating, magnetized 
object that emits a beam of electromagnetic noise (which includes radio waves, 
as well as other parts of the electromagnetic spectrum). When one looks at the 
geometry generated by a neutron star, and computes deviations from flat space of 
the metric components, one finds (on the surface of the star) 



for a star of (typical) mass 1.4Af0, and radius R — 10 km. This is a 40% 
deviation from flat space. By contrast, we recall that, in the solar system, the 
largest metric deviation from flat space occurs on the Sun's surface and is of 
order GM/{c^R) ~ 10~^. We should therefore a priori expect that such ob- 
jects might provide tests that go beyond the solar system ones in probing some 
of the strong-field aspects of relativistic gravity. In addition, in the solar sys- 
tem, the time-irreversible radiative aspects of gravity (i.e., radiation reaction ) 
are negligible (this is why we focussed above on time-symmetric interactions). 
Here, not only are strongly nonlinear effects relevant, but one must also take 
into account the time-dissymetric effects linked to using a retarded propagator 



a: 5 (t — ''^^ 1 / \xa — xb\- The corresponding time delay '^^ is typ- 



ically ~ 1 sec (since typical separations between the two objects are of order 
300, 000 km) and plays an essential role in the equations of motion of a binary 
pulsar. As a consequence, binary pulsars have given us firm experimental evi- 
dence for the reality of gravitational radiation, and for the fact that on-shell grav- 
itational radiation is described by two transverse tensorial degrees of freedom 
travelling at the velocity of light. 

In practice, only a subset of the known binary pulsars can be used for testing 
the strong nonlinear regime of GR, and/or its radiative regime. Among these, the 

''*Except in a few cases where the companion is another compact (though less compact) star rem- 
nant, a white dwarf. 



500 = -1 + 



2GM 



~ -1 + 0.4 



(3.28) 
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very best ones are PSR1913+16 (where the numbers 19hl3m and +16 deg mea- 
sure angles on the sky), PSR1534+12, PSRJl 141-6545, and PSRJ0737-3039, the 
first double binary pulsar (made of two radio pulsars, simultaneously emitting to- 
ward the Earth). 

>From the theoretical point of view, methods have been developed to deal 
with strongly self-gravitating objects, both in Einstein's theory and in alternative 
theories (see [39] for a recent review, and references). To adequately discuss the 
observations of binary pulsars, one has had to push the post-Newtonian pertur- 
bative calculation to the 2.5 PN level, i.e., to order {v/c)^. This odd power of 
the ratio u/c is Unked to the time-dissymetric, retarded nature of the propagator 
(together with some nonlinear effects). It is the first PN level where radiation 
reaction effects arise. Any experimental test of the presence of such {v /c)^ terms 
in the equations of motion is a probe of the reality of gravitational radiation. In 
addition, as mentionned above, one must carefully treat, and disentangle, the var- 
ious strong-field effects that are linked to the self-gravity of each neutron star in 
the system. 

The existing experimental tests are based on the timing of binary pulsars. Each 
time the beam of radio waves sweeps across the Earth, one observes a pulse of 
electromagnetic radiation. The data consists in recording the successive arrival 
times, say (with N = 1, 2, 3, • • • )> of these pulses. Were the pulsar fixed in 
space, these arrival times would be equally spaced in time, i.e., = to + NP, 
where P would be the (fixed) period of the pulsar. However, in a binary pulsar, 
the sequence of arrival times is a more comphcated function of the integer N 
than such a simple linear dependence. Indeed, one must take into account many 
effects: the fact that the pulsar moves on an approximately elliptical orbit, the 
deviation of this orbit from a usual Keplerian ellipse, the deviation of the or- 
bital velocity from the usual Kepler areal velocity law, the existence of various 
additional relativistic effects: gravitational redshift, second-order (relativistic) 
Doppler effect, time-delay when the electromagetnic pulse passes near the com- 
panion, radiation reaction effects in the orbital motion, etc. To compute all these 
effects, one needs to solve Einstein's equations of motion with high (~ (v/c)^) 
accuracy. 

The final result of these theoretical calculations is to derive the so-called DD 
timing formula, which gives the N*^^ pulse arrival time t at as an explicit function 
of various "Keplerian" (p^), and "post-Keplerian" (p^^) parameters, say 

tN-to = F[N;p^;p^'^]. (3.29) 

Here, the Keplerian parameters (p^) comprise parameters that would exist in a 
purely Keplerian description of the timing: the orbital period Pi,, the eccentric- 
ity of the orbit e, the time of passage at some initial periastron Tq, and some 
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corresponding angular position of the periastron ujq, and, finally, the projected 
semi-major axis x = where ai is the semi-major axis of the orbit of the 

observable puis af^ and i is the inclination angle w.rt. the plane of the sky. The 
post-Keplerian parameters (p^^) then correspond to many relativistic effects that 
go beyond a Keplerian description, namely: a dimensionless parameter k mea- 
suring the progressive advance of the periastron k — (a;)Pf,/27r, a parameter 7t 
measuring the combined second-order Doppler and gravitational redshift effects, 
possible secular variations in Keplerian parameters e, i, Pb, two parameters r, 
s measuring the "range" and the "shape" of the additional time delay that ap- 
pears when the radio waves pass near the companion, and finally a parameter Sg 
measuring the distortion of the orbit w.r.t an ellipse. By least-squares fitting the 
observed arrival times t^v'' to the above general theoretical timing formula one 
can accurately determine the numerical values of all the Keplerian parameters, 
as well as some of the post-Keplerian ones. At this stage, the determination of 
these phenomenological parameters is (in great part) independent of the choice 
of a theory of gravity. On the other hand, in any specific theory of gravity, each 
post-Keplerian parameter is predicted to be some well-defined function of the 
Keplerian parameters and of the two masses, nii and m2, of the pulsar and its 
companion. For instance, within GR the advance of the periastron is given by 

k^^{p^,m„m,) = ^^^^^, (3.30) 

where n 27r/Pb and M ~ mi + m2, while the secular variation of the orbital 
period (caused by radiation reaction effects) is given by 

(/,....,„.) .-jg ^^8-;^f-' (CM,.)V3^, <3,30 

Note that while k is proportional to 1/c^ (IPN level), the secular variation of 
the orbital period is proportional to 1/c^ (and is indeed numerically of order 
(v/c)^). The GR-predicted value for Pi, is a direct reflection of the presence of 
O ((t)/c)^) time-asymmetric radiation damping terms in the equations of motion. 
Numerically, Pb (which is dimensionless) is predicted to be of typical order of 
magnitude P^ ^ 10^^^, which seems very small, but happens to be large enough 
to be measured with good accuracy in several binary pulsars. 

The crucial point to notice is that the GR predictions (of which two are given 
here as examples) for the link between the post-Keplerian parameters and the 
masses are specific to the structure of GR, and will be replaced, in other theories 
of gravity, by different functions fc*'^''™^ (p^, mi, 1112) ^p*'^™''' (p^, 7711,7712), 



"in general, only one of the two objects, here labelled as 1, is a pulsar. 



50 



T. Damour and M. Lilley 



etc. In particular, it has been explicitly shown in various cases (and notably in 
the case of generic tensor-scalar theories where gravity is mediated both by a 
spin-2 field and a spin-0 one) that the large self-gravity of neutron stars would 
generically enter these fimctions, and drastically modify the usual prediction of 
GR, see [40]. 

To see which theory of gravity is in agreement with pulsar timing data, one 
can proceed as follows. Within each theory of gravity, the measurement of each 
post-Keplerian parameter defines a corresponding curve in the mi, TO2 plane. 
Therefore, in general, the measurement of two post-Keplerian parameters is suf- 
ficient to determine the (a priori unknown) numerical values of the two masses 
mi and m2 (as the location where the two curves intersect). Then, the measure- 
ment of any additional post-Keplerian parameter yields a clear test of the validity 
of the theory considered: the corresponding third curve should pass precisely 
through the intersection point of the first two curves. If it does not, the theory 
is invalidated by the binary pulsar data considered. By the same reasoning, the 
measurement of n different post-Keplerian parameters yields n — 2 tests of the 
underlying theory of gravity. Many such stringent tests have been obtained in 
binary pulsar observations (more precisely, nine different tests in all have been 
obtained when considering the data from four binary pulsars). Remarkably, GR 
has been found to be consistent with all these tests. Many alternative gravity 
theories have fallen by the wayside, or their parameters have been constrained 
so as to make the theory extremely close to GR in all circumstances (including 
strong-field ones). 

Let us just give two impressive examples of the beautiful agreement between 
GR and pulsar data. In the case of the original Hulse-Taylor pulsar PSR1913-H16 
the ratio between the observed value of Pf, to that predicted by GR is given by 



pobs pgal 
|j^obs /yobsj 



= 1.0026 ±0.0022, (3.32) 



where ' is a Galactic correction. The fact that this ratio is close to one corre- 



sponds to a confirmation of the relativistic force law acting on the pulsar, of the 

symbolic form f = ^1 + • • • + (^)^^ , where the crucial last term ^ {v/c)^ 

(i.e., an effect of order 10~^^) has been verified with a fractional accuracy of or- 
der 10"^. Note that this corresponds to an absolute accuracy of order 10~^^ 
compared to the leading Newtonian term GM/r"^ ! 

The timing data from the recently discovered double binary pulsar PSRJ0737- 
3039 led to the following ratio between the observed, and GR-predicted, values 
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of the post-Keplerian parameter s 



,obs 



s 

qGR fhohs 



= 0.99987 ± 0.00050. 



(3.33) 



The agreement for this parameter is at the 5x10 ^ level. 

Summarizing: binary pulsar timing data have led to accurate confirmations of 
the strong-field and radiative structure of GR. Roughly speaking, these confir- 
mations exclude any alternative theory containing long-range field^ coupled to 
bulk (hadronic) matter. 



3.3.4. Tests of gravity on very large scales 

So far, we have mainly focussed on tests of GR on spatial scales of several as- 
tronomical units (the size of the solar system), and on scales of 300, 000 km (the 
typical separation between two neutron stars). We conclude this section on ex- 
perimental tests of gravity by mentioning the existence of tests made on very 
large spatial and temporal scales. Gravitational lensing effects by galaxy clus- 
ters allow one to probe some aspects of relativistic gravity on scales ~ 100 kpc. 
Here, one is talking of the effect of the curved spacetime metric generated by 
the cluster on light emitted by very distant quasars and passing near a galaxy 
cluster containing (in addition to visible galaxies) a lot of dark matter, as well 
as some X-ray gas. Data on the temperature distribution of the X-ray gas allows 
one to directly probe the Newtonian gravitational potential U{x) of the cluster 
(without having to assume much about the (dark) matter distribution). In turn, 
the potential U{x) determines the relativistic lensing of light, via the spacetime 
metric predicted by Einstein's theory, i.e., —goo = 1 — gij = (l + ^) 6ij. 
According to Ref. [41], the agreement is of the order of 30%. This confirms the 
validity of GR on scales ~ 100 kpc. 

Primordial nucleosynthesis of light elements (e.g.. Helium, Lithium, Deu- 
terium) in the early universe depends on both the expansion rate and on the weak- 
interaction reaction rate for the conversion between neutrons and protons. Given 
that the Hubble parameter cx Gp oc GT"*, the creation of light elements 
at early times (and high temperatures T) depends on Newton's constant. The 
comparison between theoretical predictions and observations of the abundance 
of light elements typically constrains the value of G at the time of Big Bang 
nucleosynthesis, say G^^ to differ by less than O (10%) from its current value 
Qnow (-ggg Chapter 18 of the Review of Particle Physics, |tittp://pdg.lbl.gov/( l. 



°By which, one really means here fields with range larger than the distance between the two 
pulsars, i.e., ~ 300 000 km. 
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4. String-inspired phenomenology of the gravitational sector 

4.1. Overview 

>From the previous sections, one can conclude that GR is a very well confirmed 
theory so that one might be tempted to require of any future theory (and espe- 
cially string theory) that it lead to essentially no observable deviations from usual 
4-dimensional GR. For instance, one might require that all the a priori massless 
scalar fields that abound in (tree level, compactified) string theory acquire large 
masses. However, as there is yet no clear understanding of how to fit our world 
within string theory, it is phenomenologically interesting to keep an open mind 
and explore whether there exist possibilities for deviations of GR that have natu- 
rally escaped detection so far. 

String theory predicts the existence of an extended mass spectrum (g^^, (x), 
$ (x), -Bpi/, moduli fields, etc.) from which there could result some long range 
or short range modification of gravity. The existence of branes and large extra di- 
mensions could also be sources of modified gravity {e.g., KK gravity). There 
could exist short distance effects at scales of order the string scale £s which 
are observable in cosmology or in high energy astrophysics. We shall also con- 
sider possible gravitational wave signals from string-cosmology models. Finally, 
we refer the reader to the lectures by Juan Maldacena for a discussion of non- 
gaussianities in CMB data. 

A phenomenologically interesting idea (though it is not supported by precise 
theoretical arguments) is a possible breakdown of Lorentz invariance, on large 
scale physics, linked to string-scale cutoff-related effects. An example of this is 
a modified dispersion relation of the type 

E'^ = m^+f + f5i— + /32^ + --- (4.1) 
nip nip 

where mp denotes the Planck mass. One could think that because of the large 
value of the Planck mass, any such corrections to the usual dispersion relation are 
unobservable. However, there exist astrophysical phenomena, such as high en- 
ergy cosmic rays, for instance high energy 7-rays, for which such a small change 
in this relation could be observed. For example, by comparing the times-of- 
arrival of 7-rays of different energies, one has been able to place strong limits on 
the parameter Such modifications of the dispersion relation have also been 
used in the analysis of the CMB, since, in the standard inflationary model, initial 
quantum fluctuations (the seeds of today's large scale structures) arise in the deep 
ultraviolet i.e., at transplanckian scales. Note that there exist theoretical difficul- 
tiesj^with the inclusion of the (3i term (the one which is severely constrained 

"Linked to its proportionality to 1/ mp, while most theoretical models suggest a proportionality 
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BULK GRAVITY 



Figure 5. The ends of open strings are attached to a brane, giving rise to SM particles, while closed 
strings are free to propagate in the bulk. 

experimentally), while the more conventional fourth order term would be too 
small to be observed. Note that in the case of the photon, a modification on short 

scales could imply a birefringence of the vacuum as w± = |fc| ^1 ± /3;|^^ • For 
references on these issues see [42,43]. Speaking of string-inspired astrophysical 
effects, let us mention the suggestion of Ref. [44] that string theory might imply a 
violation of the usual Kerr bound on the spin of rotating black holes: J < GAP . 

Other possible predictions of string theory arise from the picture in which 
one considers the existence of branes on which (open string) SM particles are 
confined, while (closed string) gravitons are free to propagate in the bulk (FIG. 
[5]l. The extra dimensions of the bulk can then be compactified, on a Calabi- 
Yau or simply on a torus (thereby "localizing" gravity around the SM brane). 
Constraints on the size of the compactified dimensions then come either from 
the gravitational phenomenology, or from effects on SM particles. This is the 
"large" extra dimensions idea [45] which could be tested at the LHC, and so 
is of interest today. Other realizations include models with "very large" extra 
dimensions [46], but it is less clear how they are realized in string theory. In 
the Randall-Sundrum model [46], a brane can be like a defect in a bulk with a 



to 1/mp. 
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negative cosmological constant, in which case the zero mode of bulk gravitational 
waves behaves as a surface wave localized on the brane due to the discontinuity 
located at the interface of the brane with the bulk. In the DGP model [47], the 
approximate localization of bulk gravity on the SM brane is achieved through the 
interplay of two dynamics for the gravitational sector: a 5D Einstein action, plus a 
4D "induced" Einstein action, with a different value of Newton's constant, on the 
brane. Combining the two inverse propagators, the global propagator drastically 
modifies gravity on large length scales r: 



In addition, even on length scales r < L there exist modifications of usual gravity. 
Indeed, the claim is that Newton's potential is modified as [48] 



At the phenomenological level, it is interesting that (Newtonian) gravity be mod- 
ified in this way. Estimates indicate that effects are small enough to have escaped 
detection so far, but could be seen in refined solar system experiments (e.g., Lunar 
Laser Ranging). Some authors have argued that such models may have acausal 
behaviours, with, for instance, the appearance of closed timeUke curves [49]. 

Another conceptually interesting idea involves the possible existence of sev- 
eral (parallel) Randall-Sundrum branes. The confining mechanism of gravity in 
the Randall-Sundrum model is such that the wavefunction of surface gravitons 
is exponentially decaying away from the brane. If two branes are nearby, such 
quasi-confined gravitational effects can tunnel from one brane to the other via ex- 
ponentially small effects. As a consequence, the effective Lagrangian would con- 
tain two metric tensors with two gravitons, one massless, the other massive [50]. 
There are, however, theoretical difficulties with any massive gravity theory, in 
relation with the van Dam-Veltman-Zakharov discontinuity (see, e.g., [51] and 
references therein). 

4.2. Long range modifications of gravity 

It is well known that, at tree level in string theory, there exist many massless 
scalar fields with gravitational strength coupling, the so-called moduli fields. 
Phenomenologically, one would expect that having massless scalar fields at low 
energies is undesirable (Would a theory containing such massless fields not im- 
mediately fail the GR tests discussed in the section above?). General arguments 
suggest that such scalar fields should not be expected to remain massless after 




(4.2) 




(4.3) 
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supersymmetry breaking [52]. Recently, a large "industry" has been devoted to 

try to construct explicit compactification models where all moduli are fixed, and 
actually acquire very heavy masses, which is needed if inflation is to happen in 
the usual way. Here however, in the spirit of keeping an open mind, we will 
instead assume that a scalar field remains massless in the low-energy effective 
theory, and discuss ways in which it might not disagree with existing tests of 
general relativity. In other words, we suppose there exists a flat or almost flat 
direction in the total scalar potential V {(f), such that there remains a massless 
field after supersymmetry breaking. Let us mention in this respect the idea sug- 
gested, in particular, by Eliezer Rabinovici [53] that the ultimate explanation for 
the smallness of the cosmological constant might be a mechanism of spontaneous 
breaking of an underlying scale invariance. In that case, we would expect to have 
an associated massless Goldstone boson (the "dilaton", in the original sense of 
the word). 

4.2.7. The cosmological attractor mechanism 

Let us discuss here the idea of the least coupling principle, realized via a cos- 
mological attractor mechanism (see e.g., Refs [54, 55]), which can reconcile the 
existence of a massless scalar field in the low energy world with existing tests of 
GR (and with cosmological inflation). Note that, to date, it is not known whether 
this mechanism can be realized in string theory. We assume the existence of a 
massless scalar field $ (i.e., of a flat direction in the potential), with gravitational- 
strength couphng to matter. A priori, this looks phenomenologically forbidden 
but we are going to see that the cosmological attactor mechanism (CAM) tends 
to drive $ towards a value where its coupling to matter becomes naturally <C 1. 
In the string frame, we start with an effective action of the generic form 



where $ is the massless dilaton field, x is the inflaton, to which has been asso- 
ciated a simple chaotic-inflation-type potential term, with the exception that here 
is a function of In heterotic string theory for instance, Bg, B^, Bp, B^, 
and B^ are given by expansions in powers of the string coupling = e^, as 



where the first term is the tree level term, followed by an infinite series of cor- 
rection terms involving positive powers of (or non-perturbative functions of 





(4.4) 
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Qs). Switching to the Einstein frame, and redefining g^^ and the nonstandard $ 
kinetic terms according to 



Qnv 9nv = CBg ($) Qu,^, 

,2 H 1/2 



4\BgJ Bg B, 



IP = J 

(where a prime denotes d/d$), the effective action turns into 

SeS = j d^ary^ 



(4.5) 



~ 2 ~ 2 ~ 2 

I , (4-6) 



>^ (x) X^ 



+ 



with rh?p = and in which the x terms are important during inflation while 
additional terms that include the gauge fields and ordinary matter such as 



-\bp (^) Fl, -Y^jruA [Bp {xa))] 

A •' 

\J~9tiu {xa) dx^A'^x\ - Kac (<P) 



(4.7) 



are relevant in the matter dominated era. 

As we shall see, the CAM leads to some generic predictions even without 
knowing the specific structure of the various coupling functions, such as e.g., 
m^{ip),mA{BF{>fi)), ■ ■ ■ The basic assumption one has to make is that the 
string-loop corrections are such that there exists a minimum in (some of) the func- 
tions m ((/?) at some (finite or infinite) value, Lpm- During inflation, the dynamics 
is governed by a set of coupled differential equations for the scale factor, x and 
ip. In particular, the equation of motion for (p contains a term oc — j^w^ (1,5) y^. 
During inflation (i.e., when x has a large vacuum expectation value, this coupling 
drives <^ towards the special point where m-^{(fi) reaches a minimum. Once 
(fi has been so attracted near cpm, <f essentially (classically) decouples from x (so 
that inflation proceeds as if cp was not there). A similar attractor mechanism ex- 
ists during the other phases of cosmological evolution, and tends to decouple ip 
from the dominant cosmological matter. For this mechanism to efficiently decou- 
ple ip from all types of matter, one needs the special point p„i to approximately 
minimize all the important coupling functions. This can be naturally reahzed by 
assuming that (pm is a special point in field space: for instance it could be the 
fixed point of some Z2 symmetry of the T- or ^-duality type (so that one could 
say that "symmetry is attractive"). An alternative way of having such a special 
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point in field space is to assume that (p,n = +ocp^is a limiting point where all 
coupling fiinctions have finite limits. This leads to the so-called runaway dilaton 
scenario [55]. In that case the mere assumption that Bi (<i>) ~ c\~\- O (e~^*) as 
<i> +CX) implies that — +00 is an attractor where all couplings vanish. 

4.2.2. Observable consequences of the Cosmological Attractor Mechanism 
Before discussing the observational predictions of the CAM, let us remind the 
reader of a few facts that are relevant for studying the possible effects of a string- 
inspired modification of gravity. The main somxe of modification of gravity 
comes from the fact that the "moduli" field tp will influence the values of the 
masses of the (low-energy) particles and nuclei. This means that the classical 
action of, say an atom A, will be 

mA{(p)dsA = - J mA{ip)\J-gfj.udx'Xdx'X (4.8) 

where g^i^ is the Einstein-frame metric. Then, one finds that the scalar field ip 
will be coupled to the atom A with the strength oaVG, where the dimensionless 
coupling strength aA (with the same normalization as the one discussed above 
for usual tensor-scalar theorie^^ is simply given by 

= — In (93) . (4.9) 

To see better the various ways in which ip might enter into niA, let us consider 
for instance the various parts constituting the mass of an atom: 

ruA ((^) - Zmp + Nrun + Zm, + + (4. 10) 

where Z is the atomic number, the mass of the proton, iV the neutron number, 
m„ the mass of the neutron, ra^ the mass of the electron, E^^^^^ and £;™cieus 
the nuclear and Coulomb interaction energies of the nucleus, respectively. In 
addition, one must note that the mass of the proton is given by 

mp((p) = oAqcd (53 (V?)) + (V') + .4^^-) 

The main scale that determines the mass of the proton is Aqcd. It depends on all 
the moduli including the massless field (p and is roughly of the form 



Motrin.. (4.12) 



Aqcd (^) = iv) (^) exp [- 

'*This is viewed as a strong-(bare-)coupling limit, by contrast to tlie usual weak-coupling limit 
(/3 — > — 00 and 3> — > —00. 

"in particular, the effective Newton constant for a Cavendish experiment between a body made of 
atoms A and another one made of atoms B is G^^^ = G(l + a 405). 
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Here, Cg is the confonnal factor from the string to the Einstein frame. The most 

important contribution to the dependence of Aqcd is that given by the Lp de- 
pendence of the exponential term. This dependence comes from the well-known 
running (via the /3-function of SU(3)) of some (unified) gauge coupling constant 
between its value cx Bp{(f) considered at a GUT-scale cut-off (here ap- 
proximately related to Mstj-ing), to a value of order unity at the confining scale 
Aqcd- The other contributions to the mass of the proton are the quark masses, 
which are determined by the vev of the Higgs boson and by the Yukawa coupling 
constants, which, again, are expected to be functions of (p at high energy. There 
also exists a contribution from the electromagnetic sector since part of the mass 
of the proton is a function of the fine structure constant ofom if)- Finally the nu- 
clear binding energy of a nucleus is quite important and must also be expressed 
as a function of basic scales. In an approximate form it reads 

^nucleus ^^^^z)a3 + {N + Zfl^ 63 (4.13) 

where 

^, chiral limit , _2/\ 

In the chiral limit (i.e., taking the quark masses to zero) one gets a non-zero 
limit (jchirai limit ^.^ ^jjjQjj must bc added a term approximately proportional to 
the squared pion mass. In turn, is proportional to the product of Aqcd and 
tHu + TOd, both of which are expected to be functions of ^p. Incidentally, let us 
note that there exists a dehcate balance between attractive and repulsive nuclear 
interactions [56|, which implies a strong sensitivity of the binding energy of nu- 
clei to the value of the quark masses [57]. A recent result shows that if the quark 
masses were to increase by 50% (at one Itr, or 64% at 2fT), all heavy nuclei would 
fall apart because there would be no nuclear binding [58]. 

At leading order, the mass of any nucleus is a pure number times Aqcd- In 
this approximation, would depend universally on <^ (via Aqcd (</'))> and the 
scalar coupling strength cxa would be independent of the atomic species A con- 
sidered. As a consequence, there would be no violation of the universality of free 
fall. This shows that the violations of the universality of free fall will depend on 
the small fractional corrections in proportional to the ratios 

, , and aem- (4.15) 

Aqcd Aqcd 

When differentiating the mass of an atom w.r.t. ip, say 

ruA iv) = AAAqcd (l + !I^^±^ + e^ ^^^^'^" + sTa.^] , (4.16) 
V Aqcd Aqcd / 
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where J\f is a pure number (which depends on N and Z), one obtains for the 
scalar coupHng strength a a if) = ^ ™A if) an (approximate) expression of 



the form 



(4.17) 

d 

where ahad = In Aqcd if)- When the CAM has attracted near a value 

dip 

fm which minimizes all the separate coupling functions entering the various 
ingredients of mA {'p), each term in the above expression for ua {'p) will be 
(approximately) proportional to the small difference if — (pm- As a consequence 
all the contributions to ua (</?) will be small, so that all the observable deviations 
from GR will be naturally small. 

Let us describe more precisely the possible observable consequences of the 
CAM. In this mechanism, the couplings of the massless scalar field to the vari- 
ous physical sectors are not assumed to be initially small (they are given by the 
various coupling functions Bi{Lp) entering the Lagrangian, and these functions 
are "of order unity"). However, via its coupling to cosmological evolution, the 
scalar field is driven towards a point where the couplings to matter become small, 
but not exactly zero. Indeed, one can analytically estimate the "efficiency" of the 
cosmological evolution in driving ip towards ip„i, and one finds some expression 
for the differenc^^(5(/7 = ip — p^n [54,55]. The deviations from GR are all pro- 
portional to the small quantity Sip'^ because the scalar coupling strengths aA^ (Xb 
are proportional to 5ip, and all "post-Einstein" observables contain two scalar 
couplings , say a^a^ when talking about the scalar exchange between A and B 
(for instance the modified gravitational constant for a Cavendish experiment in- 
volving two bodies made of atoms A and B is Gab = G {1 + aACtB))- In addi- 
tion to predicting small values for the (approximately composition-independent) 
"post-Einstein" parameters 7 and (3 this mechanism also predicts various (small) 
violations of the equivalence principle. 

For instance, the above expressions for the ingredients entering and ua 
lead to generic predictions about the type of violation of the universality of free 
fall that one might expect in string theory. Indeed, one finds that the fractional 
difference in the free fall acceleration of two bodies (made of atoms A and B) 



"when ipm is infinite, 6ip = ip — ip,„ is replaced, e.g., by e 
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takes the form 



aA - as 



(a) 



2 X 10 



cbA 



AB 

N - Z 
M 



N ■ 



M 



AB 



AB 



with 



E 

M 



z{z-i) 

.1/3- 



(4.18) 



(4.19) 



{N + zy 

where (AQ)^^ = Qa — Qb< and where the first, second and third terms in 
the brackets are contributions from the Coulomb energy of the nucleus (oc Ofem), 
and from the (^-dependence of the sum and difference of the quarks masses, i.e., 
rriu + nid and m„ - rrid. 

This mechanism also predicts (approximately composition-independent) val- 
ues for the post-Einstein parameters 7 and parametrizing IPN-level deviations 
from GR. They are of the form 



7 =-2: 



a 



had 



~ -2a? 



had' 



''had 



and 



1 «: 



had 0<p 



1 



da 



had 



In this model, one in fact violates all tests of GR. However, all these violations 
are correlated. For instance, using the numerical value A (^) ~ 2.6 (which 
applies both to the pair Cu-Be and to the pair Pt-Ti), one finds the following link 
between equivalence-principle violations and solar-system deviations 



(4.20) 



(4.21) 



Aa 



-2.6 X 10-^7. 



(4.22) 



Given that present tests of the equivalence principle place a Umit on the ratio 
Aa/a of the order of 10~^^, one finds I7I < 4 x 10^*. Note that the upper hmit 
given on 7 by the Cassini experiment was 10~^, so that in this case the necessary 
sensitivity has not yet been reached to test the CAM. 

As another example, one can compute the evolution of the fine structure con- 
stant w.r.t. time. Given that it is a function of (f, and that (p evolves as a function 
of cosmological evolution due to its coupling to matter, a^ni is indeed a function 
of time, and its time derivative can be written as 



37 in "em 

at 



±10 



-16 



l + «0 



30^ 
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, Aa 



(4.23) 



String theory, gravity and experiment 



61 



The first square root on the r.h.s. of this equation can also be written as 

^mO^m H~ 4ri^ct^ 24) 

where ilm denotes the fraction of the cosmological closure density due to dark 
matter, and the scalar coupling to dark matter, while il^ and Uy denote the 
corresponding quantities for "dark energy" (or "vacuum energy"). For instance, 
if we assume ^ 1 (so that (p is a kind of "quintessence") while am ^ 1, we 
see from the result above that the current experimental limit ^ < 10^ im- 
plies the following upper bound on a possible time variation of the fine-structure 
constant: (Inacm) < 10"^^ yr^^- This upper bound is below the current 
laboratory limits on a /a, but comparable to the Oklo limit mentioned above. 
When working out the generic predictions of the runaway dilaton version of the 
cosmological attractor mechanism, one finds that it naturally predicts (when as- 
suming an inflationary potential cx x^) a level of deviation from GR of order 
—7 ~ 4C X 10^^, corresponding, for instance, to a violation of the equiva- 
lence principle at the level Aa/a ^ C x 10^^^. Here, C is a combination of 
model-dependent dimensionless parameters, which are generically expected to 
be "of order unity". This suggests (if C is smaller, but not much smaller than 
1) that the current sensitivity of equivalence principle experiments may be close 
to what is needed to test the deviations from GR predicted by such a runaway 
dilaton. Let us note in this respect that ongoing improved lunar laser ranging ex- 
periments will probe Aa/a to better than the ^ 10^^"^ level, and that the CNES 
satellite mission MICROSCOPE (to be launched in the coming years) will reach 
Aa/a ^ 10~^^. Another more ambitious satellite mission (which is not yet 
approved), STEP (Satellite Test of the Equivalence Principle), plans to probe 
violations of the equivalence principle down to the 10~^* leve]^ In addition, 
post-Newtonian solar system experiments at the 10~^ level would be of interest. 
The approved micro-arcsecond global astrometry experiment GAIA will probe 
7 ~ 10~^, while the planned laser experiment LATOR might reach 7 ~ 10~^. 
In addition, the comparison of cold-atom clocks might soon reach the interesting 
level ^ (In aenO'-lO-i^ yr-^. 

Finally, let us mention that one can combine the basic mechanism of the CAM 
(which consists in using the coupling of (p to matter, i.e., the presence of a term 
of the form a((^)/9mattor in the action) with the presence of a "quintessence"- 
like potential V{(p) oc This yields the "chameleon" mechanism [60] in 

which both the value towards which ip is attracted, and the effective mass 
(or inverse range) of (p, depend on the local matter density Pmattor- Whatever be 

^'Let us also mention the suggestion [59] that atom interferometry might be used for testing the 
equivalence principle down to the Aa/a ~ 10"^'' level. 
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one's opinion concerning the a priori plausibility of having some nearly massless 
moduli field surviving in the low-energy physics of string theory, it is clear that 
such experiments are important and could teach us something new about reality]^ 



5. String-related signals in cosmology 

5.7. Alternatives to slow-roll inflation 

In the usual inflationary scenario, the period of exponential expansion is based 
on the slow roll mechanism, i.e., one has to assume a sufficiently flat potential 
so that the scalar field, the inflaton, slowly rolls down to its minimum in such a 
way that the approximate equality ~ —p^p lasts sufficiently long, say for a 
minimum of 60-70 e-folds. The simplest inflationary Lagrangian reads 

C = -^{d^f-V{^) (5.1) 

with a usual kinetic term and a potential V {ip). This simple inflationary frame- 
work leads to specific predictions such as a relation between the ratio of tensor 
to scalar primordial perturbations and the "distance" in field space over which 
(p runs during inflation (the so-called Lyth bound, see the lectures by Juan Mal- 
dacena in these proceedings). Let us, however, emphasize that these predictions 
(which lead to constraints on the model) do depend on the assumption that infla- 



tion is realized by the simple action (5.1 1 with a slow -roll potential. There are. 



however, other ways of realizing inflation, in which these constraints might be re- 
laxed. Let us note in this respect that inflation can be realized even if the potential 



V (if) in (5.1 1 is not of the slow-roll type [61]. Moreover, one may have inflation 
without a potential at all if the Lagrangian is a complicated enough function of 
X = — {dipf . Indeed, if one has an action of the type C = p{X), one finds 
that there can exist attractors toward a de Sitter expansion phase, corresponding 
to a line where the effective equation of state deduced from C — p{X) isp = —p 
{e.g., k-inflation [62]; ghost inflation [63]). To have a "graceful exit" from this de 
Sitter phase one needs, for instance, to introduce some additional tf dependence 
in £. It has been suggested in [64] that such a mechanism might be realized in 
string theory, via a Dirac-Born-Infeld-type action, say 

piX,^) = -^^f^-iyVi^). (5.2) 

^-For instance, if one considers it very unlikely that such a field can exist, these experiments are 
important because they can falsify string theory. By contrast, if one finds a violation of the equivalence 



principle which is nicely consistent with the prediction l |4.18| for the composition dependence of a 
moduli field, this might be viewed as a confirmation of string theory. 
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In such a "DBI inflation", the use of non-standard kinetic terms greatly relaxes 
the restrictions imposed on the flatness of the potential V {(p) which must be im- 



posed in the usual case of (5.1 1. It also tends to produce larger non-gaussianities 
in the CMB [65]. Let us also point out that the use of a non-linear kinetic term 
might significantly affect the Lyth bound. For instance, if one considers the ac- 
tion 

C^p{X,^)^K{X)~V{ip), (5.3) 

where K {X) is a non-linear function of the kinetic term X = — {dtpf, one finds 
the following modified form of the relation between the ratio r of tensor to scalar 
primordial perturbations and the derivative of tp w.rt. the number of efolds N: 



(5.4) 



Here a is an additional amplification factor, which is given by the following ex- 
pression in terms of the kinetic function K {X) 

, lforii:=— , 
2K' 2 ' 

= < 1 for DBI type : - Vl -X, 



^' ' > 1 for, e.g., - — (1 - X)" , with a<-. 

2a 2 



K 



A large amplification factor a ^ 1 would (formally) correspond to a relaxed 
Lyth bound on the excursion of Lp, given a minimum number N of efolds. It is 
interesting to note that a — \ (unchanged Lyth bound) in the DBI-type square 
root model. However, we note that a more general power a, with a < 1/2, would 
formally relax the Lyth bound. [A more detailed study is, however, necessary for 
seeing whether the bound is physically relaxed.] 

The present section presented only a very partial and sketchy picture. It was 
only intended as an illustration that folklore results and constraints on inflationary 



models do depend on using the standard slow -roll action ( 5. 1 1, and that there exist 
other mechanisms in which those results and constraints might be different and 
possibly relaxed. 



5.2. Cosmic superstrings 
5.2.1. Phenomenological origin 

The existence and detection of cosmic superstrings is an exciting possibility that 
was first suggested in Ref. [66], then kept alive for a number of years, and re- 
cently revived dramatically notably in Ref. [67] and in other papers [68-70]. 
They arise in brane antibrane scenarios where the inflaton is the brane-antibrane 
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Figure 6. Left: The brane-antibrane distance as a scalar field ip. Right: V (ip) behaves as ci — C2ip ^ 
for large brane-antibrane separations 



separation (see e.g., [71]). In such scenarios, for large enough brane-antibrane 
separations, the potential behaves like ci — cij^'^ such that it satisfies the slow 
roll conditions (F1G.[6|. When the branes are near, some of the modes connect- 
ing the two branes become tachyonic, i.e., a complex field (with kinetic term 
—dTdT) having a potential V^(|Tp) with wrong-sign curvature near T = 0. 
This instability can generate topological defects since the phase of the vev of 
T need not be uniformly the same all over space. Contrary to the situation in 
which strings are created at the beginning of inflation and then diluted away, this 
scenario naturally produces strings at the end of inflation so that they are not di- 
luted by the expansion. For causality reasons, the value of the field's vev in a 
given Hubble patch should be uncorrected to that in other Hubble patches. This 
is what creates a network of strings and one can then compute the initial den- 
sity and correlation length of the string network. The string tension ji in Planck 
units, i.e., the dimensionless parameter G ji, where G is Newton's constant, was 
initially thought to be high, of the order of lO^'^ at best, because of the string 
theoretic origin of these objects and of the then expected relation between a' and 
the Planck length. However, in models with warping factors and large fluxes, 
the string tension can be lowered to much smaller values. In practice, the string 
tension is tuned to fit current CMB data. Tye and collaborators [68] find a win- 
dow of the type 10"^^ < G /i < 10"'^, while in die more detailed KKLMMT 
model [67] one finds G^i - 10"^". 

In trying to gain insight into the observational predictions that can be made 
from cosmic superstring models, one must consider not only the stretching by the 
cosmological expansion of an initial network of cosmic strings with a correlation 
length of the order of the Hubble scale , but also string interactions. A string can 
for instance self-intersect or two strings can intersect and reconnect. The Hubble 
expansion tends to locally straighten out the strings while interconnections tend 
to produce loops and small-scale structure. Given an initial correlation length and 
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reconnection probability p, working out the time evolution of a string network is 

essentially a classical problem. Two types of strings develop, long strings with 
correlation length of the order of the time scale t, and small loops that loose 
energy by gravitational radiation [72-74]. 

In order to define the typical size of string loops (at the time they are formed) 
we introduce a dimensionless parameter a such that ^loop {t) = etc t. It was 
initially thought that a = 50 G /U, an estimate Unked to the idea that gravita- 
tional damping is the essential mechanism wich determines the lifetime of loops. 
More recently, is has been suggested that a might be significantly different from 
50 G/U. There is, however, no consensus on the "correct" value of q. Estimates 
vary between a ^ (50 G fi)^^, with /? > 1 (leading to "small loops") and a 0.1 
(leading to large loops). [For an introduction to this problem, and references, 
see e.g., the talk of Joe Polchinski at the 2007 String meeting in Madrid.] Hap- 
pily, some of the predictions we shall discuss below (notably those concerning 
the observabihty of gravitational waves from a cosmic string network) are rather 
insensitive to the value of a. 

Several numerical simulations confirm the tendency of string networks to dis- 
play a scale-invariant behavior [75, 76] . There have been recent attempts at refin- 
ing the theoretical description of string networks [77,78]. However, there is, to 
date, no consensus among experts as to the typical size distribution of loops (i.e., 
the dominant gravitational wave-emitting string type). In several simulations, the 
distribution of the size of loops is bimodal, with one peak at a ^ 0.1 and another 
peak at the UV cutoff. It has been argued by Vilenkin and collaborators that only 
the "large loop" part, i.e., a ~ 0.1, will survive. 

In the following, assuming KKLMMT-type brane inflation and the stability of 
strings over cosmologically interesting time scales, we discuss the phenomeno- 
logical predictions made by treating p and a as free parameters and their possible 
observable signals. 

5.2.2. Observational signatures 

Partly for historical reasons, the phenomenology of cosmic strings has been stud- 
ied mostly in the context of CMB observations. Slow-roll inflation generates a 
random ^ angular distribution on the sky that fits well the observations. Adding 
a random network of cosmic strings generates additional (non-Gaussian) fluctua- 
tions in the CMB which have less angular structure (the string has a lensing effect 
proportional to its velocity v over the sky, 5T/T ^ SttGuv^). CMB observa- 
tions can then be used to place an upper bound on G/i, of order 5 x 10~^. Much 
smaller values of the string tension will not lead to any observable signature in 
the CMB. Let us also mention that cosmic (super)strings might be detected via 
their gravitational lensing of galaxies, or microlensing of stars. 

By contrast to the CMB (or lensing) observations, which are only sensitive to 
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String tensions Gii > 10^^, existing or planned gravitational wave interferom- 
eters could detect cosmic (super)strings with tensions in the much wider range 
10~^^ < G fi < 10^^. Let us recall that a gravitational wave (GW) detector is 
actually measuring tidal forces, and more precisely a component of the Riemann 
tensor projected "along" a detector having a quadrupolar structure]^ In other 
words, a GW antenna measures the second time derivative h (t) of a projection 
of the metric fluctuation /i^j,. Current detectors are sensitive down to the level 
h ~ 10~^^, for frequencies / ~ 100 Hz. 

The GW signal from a string network is an incoherent background of GWs 
made of the superposition of all GWs ever emitted by string loops (from zero 
to very large redshifts). This signal is distributed over a very large spectrum 
of frequencies (including wavelengths of the size of the universe, as well as very 
short ones). In order to determine the frequency distribution, the number of loops 
and how they evolve, one needs to know the evolution of the universe during the 
inflationary, radiation and matter dominated eras. 

Besides detecting the GWs from a string network in a man-made interferom- 
eter, another observational possibility lies in the timing of isolated millisecond 
pulsars. In a stationary spacetime, the pulses emitted by an isolated pulsar would 
be observed on Earth (after correcting for the Earth motion) at very regular in- 
tervals. By contrast, in presence of a fluctuating background of GWs, the times 
of arrival of successive pulses would fluctuate, and exhibit some red noise. Pul- 
sar timing over some time interval T (which is typically several years) is most 
sensitive to the part of the GW frequency spectrum with frequencies / ~ 1/T. 
Therefore, pulsar timing is most likely to detect long wavelength GWs (several 
light years long). 

Along with LIGO-type ground based interferometers, a space-based one, the 
Laser Interferometer Space Antenna (LISA) has been conceived, with arm lengths 
of the order of 10^ km instead of the 3 or 4 km ones constructed on the ground. 
LISA can therefore explore much smaller frequencies. The best achievable sen- 
sitivity for LIGO-type instruments is reached for frequencies / ~ 100 Hz, i.e., 
rather fast events lasting ~ 10^^ seconds, while space experiments may probe 
events with periods ^ 1000 sees, which are quite slow events (see FIG. [7]). In 
Ref. [74], the possible existence of sharp gravitational wave bursts above the 
background caused by string cusps was pointed out. For a typical oscillating 
loop, there occurs a cusp once or twice per oscillation with the extremity of the 
cusp going at the velocity of light and emitting a strong gravitational wave sig- 
nal in the direction in which the cusp is moving. Statistically, these events are 
random. A GW burst will be detected if it happens to be emitted towards the 



This is tlie spin-2 analog of saying that electromagnetic antennas are sensitive to the projection 
of the electric field along the direction of a dipolar antenna. 
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detector. Under some conditions, these cusps can create signals which stand 
much above the quasi-Gaussian random mean square background "GW noise". 
This raises the exciting possibility that LIGO/VIRGO/GEO or LISA might detect 
GW signals emitted by giant superstrings at cosmological distances. 

Let us now give an introduction to the physics behind the occurrence of those 
cusps, and the associated emission of GW bursts. 

5.2.3. String dynamics 

We consider the string position as a function of the worldsheet coordinates r 
and (7. We treat the string dynamics in a locally flat spacetime. Introducing the 
lightcone coordinates in conformal gauge, 

a± = T ± a, (5.5) 

X'^ (r, a) satisfies 

i-^i-_X^M=0 (5.6) 
such that the generic string solution is the sum of left and right movers 

X^{T,a) = ^[X't{<j+) + X^{a^)] (5.7) 

in which the factor 1/2 is introduced for convenience. The Virasoro constraints 
read 

(d+xt)' = 

= ^^-^^ 
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In the time gauge, the worldsheet is sHced by constant coordinate time hyper- 
planes = x*^ = r, so that 

X°(T,a) = r= + (5.9) 

We thus have = (7+ and X^ = a". Then d±X'^ = 1 contributes a —1 in the 
Virasoro constraints, so that 

(d+X^Y ^ -1 + {d+X')\ (5.10) 

and similarly for the — equation. This means that the derivatives (w.r.t. their 
argument) of the spatial components Xj.{a±) of the left and right modes are 
constrained to be unit euclidean vectors: 



' = 1. (5.11) 



The Xj. are periodic and the time derivative of the spatial component of X are 
unit vectors. We may now use a representation first introduced in Ref. [79]. The 
derivatives X^. can be seen as drawing two curves on the unit sphere (the Turok- 
Kibble sphere). In addition, as X^. is periodic in three-dimensional space (there 

is no winding), we have J da±Xj. = 0. As a result, the "center of mass" of both 
left and right moving curves must be at the center of the sphere. This implies 
that the two curves genericall^p^ intersect twice [80]. Now, the main point is 
that an intersection between the two curves represents a cusp. More technically, 
such an intersection corresp onds to particular points on the string worldsheet 



at which the two null (see (5.8 1) tangent vectors Xf^ and X^ are parallel in 



spacetime. In general, the string worldsheet intersects locally the light cone along 
the two distinct directions X'^ and X^ . The cusps are special points where the 
worldsheet is tangent to the lightcone (see FIG. [8]l. This is a singularity of the 
classical worldsheet at which a strong gravitational wave signal is emitted along 
the common null vector. Let us now indicate how one computes the emission of 
gravitational wave bursts from cuspy strings. We consider Einstein's theory in 
the linearized approximation, 

gt^u{x) ^ri^i, + h^,y{x) . (5.12) 

We use the harmonic gauge, d'^h^i, = 0, so that Einstein's equations simplify to 

□ V = -le^GT^. (2^) , (5-13) 



There exist, however, specially contrived curves that can avoid intersecting. 
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Figure 8. Left: The lightcone generically intersects tlie worldslieet into two separate null directions 
corresponding to the velocity of the left and right movers. Right: When a cusp occurs, the two null 
vectors are parallel and the worldsheet is tangent to the lightcone. There results a large burst of 
outgoing gravitational radiation. 



where 



1, 



(5.14) 



The stress-energy tensor is obtained by differentiating the Nambu action w.r.t. 
f/pi,. Taking its Fourier transform, one finds 



T'^" (fc^) 



lk.{X++X^) 



(5.15) 



where {^v) indicates symmetrization over the indices /i, and where, in the 
exponential, we have replaced X^^ by the half-sum of the left and right movers. 
The fundamental period of a loop of length £ is Ti = ^ . Note that £ is the invariant 
total length where n is the string tension. In string theory, one usually uses a 
worldsheet gauge where £ is either 1 or 27r but here one finds it more convenient 
to use a gauge where a and r are connected to an external definition of time 
(namely X^ = x*^ = r). 

We wish to compute the integral giving T^" (fc^) over a periodic domain 
in the r , a plane. We can rewrite the integral as an integral over d(T+d(T_. This 
yields the famous left-right factorization of closed string amplitude^ and the 
Fourier transform of the string stress-energy tensor reads 



(5.16) 



^""Though we are doing here a classical calculation, one recognizes that the result is given by the 
graviton vertex operator. 
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where 

I±= I dcr±X^{a±)e-i''-^^. (5.17) 



By solving Einstein's equation (5.13 i, one finds that the spacetime-Fourier trans- 
form of the source on the r.h.s. actually gives the time-Fourier transform of the 
asymptotic GW amplitude (emitted in the direction = k'^/k^), i.e., the time- 
Fourier transform of the quantity {t — r, ft) appearing in the asymptotic ex- 
pansion 

Here the 1 /r decrease in amplitude as a function of distance away from the string 
is caused by the retarded Green's function in 3+1 dimensions, k^j^ is a function of 
both the time variable and the angle of emission. As we just said, the time -Fourier 
transform of k^'^ is proportional to the spacetime-Fourier transform T^"^ {k) of 
the source, and is explicitly given by 

= |/|y'dte2-/(*-'-)«:^^(t-r,n) ^^^^^ 

= 2G^JL\f\I'f {uj,ijjn)l'l\uj,ujn) . 

This formula shows that we can compute what is observed in a GW detector as a 
function of string tension, frequency, and the product of two integrals involving 
left and right moving modes. 

We can then estimate the generic features of the GW burst emitted by a cusp 
by noticing that, in the Fourier domain, each integral I'^ is dominated (when 
considering large frequencies: / ^ T^^) by the singular behaviour of the two 
integrands X'^ {cr±) e^^'^-^^ near a cusp. The calculation proceeds by (Taylor) 
expanding the vectors X!^ and X!^ in powers of a± . One finds that the first few 
leading terms in this expansion can be gauged away, so that the signal amplitude 
is much smaller than what could have (and had) been initially thought. After 
Fourier transforming back to the time domain, it is finally found that [74] 



k (t) CX \t ~ tf. 



5/3' (5-20) 



As this result seems to crucially depend on the presence of a mathematically sin- 
gular behaviour of the classical string worldsheet at a cusp, one might worry that 
quantum effects could blur away the sharp cusp, and make the above classical 
burst signal disappear It was checked that this is not the case [81] (the basic rea- 
son being that, finally, the strong GW burst signal is emitted by a large segment 
of the string around the cusp). 
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5.2.4. Gravitational waves from a cosmological string network 
In order to understand the observational signature of a cosmic string network and 
not just a single string, one must combine the analysis of the previous section 
with the cosmological expansion of a Friedman-Lemaitre universe and with an 
integration over redshift. A crucial point is then to estimate the number density of 
string loops. This density can be analytically estimated as a function of the string 
parameters, such as the reconnection probability p, and the string tension Gjj,. 
Note that the reconnection probability is expected to be quite different for cosmic 
supertrings compared to the traditionally considered field-theory strings. Field 
theory strings are expected to reconnect, when they cross, with essentially unit 
probability (p ~ 1), while fundamental or D-strings are expected to reconnect 
with a smallish probability, 10^"^ < p < 1 [82] (because of the presence of the 
string coupling, and other factors). 

The loop number density can be approximately estimated as [83] 

where the first term on the r.h.s. comes from loops that were created at redshifts 
< 1, while the ' . . . ' denotes a possible additional contribution from high-redshift 
strings. When the loop-size parameter a is smaller or equal to the "tradition- 
ally expected" value SOG/i, the contribution from high-redshift strings is neg- 
ligible (because strings decay in less than a Hubble time). By contrast, when 
a 3> SOG/i, the strings survive over many Hubble times, and the contribution 
of high-redshift strings starts to dominate the loop density. Note the somewhat 
unexpected feature displayed by the first term in n^, namely that it increases both 
as G /i and/or p are decreased. This feature is one of the features which allow 
GW signals from strings to be detected down to very small values of the string 
tension (contrary to CMB effects). Indeed, as G /i is decreased, though each 
individual string signal will decrease proportionally to G fi, there will be more 
emitting loops. After integrating over redhifts, one finds that the observable sig- 
nal is a complicated, non monotonic function of G fi. The numerical estimates 
of Ref. [83] considered the case in which the loop size parameter a < 50G/X, 
in which case the first term in ( |5.21[ ) is dominant. If, on the other hand, one as- 
sumes a ^ 0.1 (as is suggested by some numerical simulations [75]) , strings 
survive longer so that higher redshift contributions are non-negligible. It has 
been found that in such cases these contributions increase the number of loops 
(which increases the GW signal) but tend to drown the cusp signal within the 
quasi-Gaussiam random-mean-square GW background [84]. 

Based on current detector capabilities and on the sensitivity estimates for fu- 
ture detectors, one finds that if a < 50G/i, LIGO could detect G/x > 10^^^ 
while LISA could detect G fi > IQ-^'^. On the other hand, if a > 50G^ LISA 
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could reach G n> 10~^^. One has looked in the current LIGO data for the pos- 
sible presence of a background of GW's, but without success so far [85]. The 
best current bound on G /z comes from pulsar timing [86] and is roughly at the 
G n < 10 ~^ level (which is about three orders of magnitude more stringent than 
the limits than can be obtained from CMB data). 

Gravitational wave detectors are thus excellent probes of cosmic (super)strings. 
There is therefore the possibiUty that they could confirm or refute KKLMMT- 
type scenarios in a large domain of parameter space. However, there are large 
uncertainties in string network dynamics which prevent one from being able to 
make reliable analytical estimates. If one is in a region of parameter space where 
the rather specific cusp-related signals are well above the r.m.s. background one 
might find rather direct experimental evidence for the existence of cosmic strings. 
There would however remain the task of discriminating between string theoretic 
strings and field theoretic ones. One way would be (assuming one could strongly 
reduce the string network uncertainties) to determine the reconnection probability 
p from its influence on the loop number density, and, thereby, on the recurrence 
rate of observed signals. Another more ambitious possibility would be to exploit 
the presence of two populations of strings, namely D and F strings, in D-brane 
anti D-brane annihilation, and attempt to measure two different values of G /U, 
the ratio of which satisfies = lipids- 

6. Conclusion 

We hope that these lectures have shown that gravity phenomenology is a poten- 
tially interesting arena for eventually confronting string theory to reality. 
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